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Abstract 

We study the two-loop dilatation operator in the noncompact SL{2) sector of QCD and su- 
persymmetric Yang-Mills theories with M = 1,2,4 supercharges. The analysis is performed for 
Wilson operators built from three quark/gaugino fields of the same helicity belonging to the fun- 
damental/adjoint representation of the SU(3)/ SU(N C ) gauge group and involving an arbitrary 
number of covariant derivatives projected onto the light-cone. To one-loop order, the dilatation 
operator inherits the conformal symmetry of the classical theory and is given in the multi-color 
limit by a local Hamiltonian of the Heisenberg magnet with the spin operators being generators 
of the collinear subgroup of full (super) conformal group. Starting from two loops, the dilatation 
operator depends on the representation of the gauge group and, in addition, receives corrections 
stemming from the violation of the conformal symmetry. We compute its eigenspectrum and 
demonstrate that to two-loop order integrability survives the conformal symmetry breaking in 
the aforementioned gauge theories, but it is violated in QCD by the contribution of nonplanar 
diagrams. In SYM theories with extended supersymmetry, the M— dependence of the two-loop 
dilatation operator can be factorized (modulo an additive normalization constant) into a multi- 
plicative c-number. This property makes the eigenspectrum of the two-loop dilatation operator 
alike in all gauge theories including the maximally supersymmetric theory. Our analysis suggests 
that integrability is only tied to the planar limit and it is sensitive neither to conformal symmetry 
nor supersymmetry. 
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1 Introduction 



There is a growing amount of evidence that four- dimensional gauge theories possess a hidden 
symmetry. The latter is not seen at the level of the classical Lagrangian and manifests itself 
through remarkable integrability properties of effective dynamics of the quantum theory. The 
symmetry has been first discovered in studies of high-energy asymptotics of scattering amplitudes 
in multi-color QCD PJ 0] which led to the identification of the evolution equation for the partial 
waves with the Schrodinger equation for the celebrated Heisenberg magnet with spins being 
generators of the SL(2, C) group. Later, similar integrable structures have been observed in 
the renormalization group (Callan-Symanzik) equations for multi-particle distributions in QCD 
I3HUE]- These distributions are given by (nonperturbative) matrix elements of composite Wilson 
operators built from elementary quark and gluon fields and an arbitrary number of covariant 
derivatives in the totally symmetric representations of the conformal group 5*0 (4, 2). The scale 
dependence of the distribution amplitudes is driven by the anomalous dimensions of these Wilson 
operators which, in their turn, are defined as eigenvalues of the dilatation operator in QCD. It 
was found that the one-loop dilatation operator in QCD in the sector of the so-called maximal- 
helicity operators can be mapped in the multi-color limit into a Hamiltonian of the Heisenberg 
magnet with the spin operators in the chain sites defined by the generators of the 'collinear' 
SL(2, R) subgroup of the full conformal group jSHHE]- The phenomenon of integrability is not 
merely of academic interest as it has a number of important phenomenological applications in 
QCD (for a review see P). 

The above mentioned integrable structures are not specific to QCD. They are also present 
in generic gauge theories including supersymmetric Yang-Mills models with M = 1,2,4 super- 
charges. Supersymmetry enhances the phenomenon by extending integrability to a larger class 
of observables. In this context, the maximally supersymmetric Af = 4 Yang-Mills theory is of 
a special interest with regards to the AdS/CFT correspondence [Zj. The gauge/string duality 
establishes the correspondence between this superconformal Yang-Mills theory and the type IIB 
string theory on the AdSsxS 5 background 7\. In particular, the anomalous dimensions of com- 
posite Wilson operators in M = 4 gauge theory are mapped into the energies of string excitations 
in the dual string picture [HI 13 HH UJ] • The string sigma models on the anti-de Sitter space are 
known to possess an infinite set of conserved currents and thus turn out to be completely inte- 
grable. This was shown classically to be the case for the critical AdSsxS 5 string action [T2*l IT3] 
as well as for noncritical sigma models with the anti-de Sitter as a factor of the target space 
|14j . which are dual to particular (non-) supersymmetric Yang-Mills theories. The gauge/string 
duality hints that these structures should manifest themselves through hidden symmetries of 
the spectra of anomalous dimensions of Wilson operators in dual gauge theories to all orders of 
perturbation theory. 

At present, all studies firmly point out that the one- loop dilatation operator in certain sectors 
of QCD and supersymmetric Yang-Mills (SYM) theories with < 4 supercharges possesses 
integrability properties. In maximally supersymmetric theory the phenomenon gets extended to 
all Wilson operators . A natural question arises whether integrability persists to higher loops 
as suggested by the gauge/string duality. Recently a great deal of activity has been devoted 
to the test of higher loop integrability in the M = 4 theory and significant evidence has been 
gathered supporting its presence in various closed 'compact' subsectors specific to this particular 
theory (see Ref. JE] for a review and Ref. J7j for recent developments). These sectors embrace 
composite operators with a nontrivial R— charge for which the corresponding mixing matrices 
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have a finite size determined by the number of constituent fields which build up the former. 
They do not comprise however Wilson operators involving an arbitrary number of covariant 
derivatives. 1 For such operators, the mixing matrices could have an arbitrary large size and 
they define the 'noncompact' sector. In distinction with the compact sectors, which are intrinsic 
to models with a nontrivial R— symmetry group, the noncompact ones are present in all gauge 
theories including nonsupersymmetric ones. This offers one an opportunity to compare the 
properties of the dilatation operators in various gauge theories in the same noncompact sector 
and identify the origin of the symmetry responsible for the existence of integrability. 

In this paper we will study the dilatation operator in the noncompact SL(2) sector for the 
aligned-helicity fermionic operators in QCD and supersymmetric Yang-Mills theories with M = 
1,2,4 supercharges. 2 These operators are constructed from quark/gaugino fields of helicity +~ 
belonging to the fundamental /adjoint representation of the gauge group and an arbitrary number 
of covariant derivatives projected onto the light-cone. To one-loop order, both in QCD and all 
SYM models in the multi-color limit, the mixing matrix for such operators can be mapped into 
a Hamiltonian of the Heisenberg spin chain. The length of the spin chain equals the number 
of fermions which form a given Wilson operator and the spin operators are determined by the 
generators of the collinear SL(2; R) subgroup of full (super) conformal group. 

It is well-known that the conformal symmetry is broken in QCD and SYM theories with 
M < 4 supercharges while in the maximally supersymmetric M = 4 model it survives on the 
quantum level. However the conformal anomaly modifies anomalous dimensions starting from 
two loops only and, therefore, the one-loop dilatation operator inherits the conformal symmetry 
of the classical theory. Starting from two-loop order, the dilatation operator in the SL(2) sector 
acquires several new features. First, it receives conformal symmetry breaking corrections arising 
both due to a nonzero beta-function and a subtle symmetry-violating effect induced by the 
regularization procedure. 3 Second, the form of the dilatation operator starts to depend on the 
representation of the fermion fields, i.e., fundamental SU(3) in QCD and adjoint SU(N C ) in 
SYM theories. The difference between the two is that it is only in the latter case that one 
can select planar diagrams by going over to the multi-color limit, while in the former case the 
large- N c counting is inapplicable and the two-loop dilatation operator receives equally important 
contributions from both planar and nonplanar Feynman graphs. Thus, by studying the two-loop 
dilatation operator in the SL{2) sector we are able to identify what intrinsic properties of gauge 
theories are responsible for the existence of the integrability phenomenon per se: 

• Conformal symmetry 

• Super symmetry 

• Planar limit 

We are able to give definite answers to all of the raised questions by performing the two-loop cal- 
culation of the anomalous dimensions of the aforementioned aligned-helicity fermionic operators 
in QCD and all SYM theories. The computation is performed within two different renormal- 
ization schemes based on dimensional regularization and dimensional reduction. The latter is 

lr The exception is the study in Ref. |18| of spin-three scalar operators with three covariant derivatives. 
A preliminary account of this work was published in a letter [T5|. 
3 It worth mentioning that in distinction with the compact sectors, the dilatation operator of gauge theories in 
the SL(2) sector defines one of the component of the spin operator on the spin chain side. 
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required for preserving supersymmetry in SYM theories. Actually, we find that the anomalous 
dimensions in the two schemes only differ by an additive constant indicating that the integrability 
phenomenon (if present) is scheme independent. 

Our finding can be summarized as follows. According to Liouville, the dynamical model 
is integrable provided it acquires enough integrals of motion to match the number of degrees 
of freedom. For an all-loop dilatation operator H(A), depending on 't Hooft coupling constant 
A and acting on a Wilson operator built from L constituent fields and an arbitrary number 
of covariant derivatives, integrability would require, in general, the existence of L conserved 
charges. Two of the charges — the light-cone component of the total momentum of L fields and 
the scaling dimension of the operator — follow immediately from Lorentz covariance of the gauge 
theory. 4 However, the identification of the remaining charges <2%(A) with k = 3, . . . , L is an 
extremely nontrivial task. The eigenvalues of the charges qt define the complete set of quantum 
numbers parameterizing the eigenspectrum of the dilatation operator. Integrability imposes a 
nontrivial analytical structure of anomalous dimensions of Wilson operators and implies the 
double degeneracy of eigenvalues with the opposite parity 1201 HI]- To the lowest order of 
perturbation theory this property has been verified in Refs. [3J H] with the help of the Quantum 
Inverse Scattering Method [21]. In particular, it was found there that breaking of integrability 
leads to lifting of the degeneracy in the eigenspectrum of the one-loop dilatation operator. 5 We 
argue that the same relation between integrability and degeneracy of the eigenstates holds true 
to two loops and allows one to identify high loop corrections to the conserved charges. 

Going over through diagonalization of the two-loop dilatation operator, we find that the 
desired pairing of eigenvalues occurs for three-gaugino operators in SYM theories with M = 1,2,4 
supercharges and the SU(N C ) gauge group. In that case, nonplanar diagrams provide vanishing 
contribution and, as a consequence, the exact anomalous dimensions coincide with those in the 
multi-color limit. The three-quark (baryonic) operators in QCD are well-defined for the SU(3) 
gauge group only. Due to different color representation, the two-loop dilatation operator in 
QCD receives corrections from certain nonplanar diagrams which yield breaking of pairing of 
eigenvalues. To trace the origin of breaking of integrability we split the expression for the two- 
loop anomalous dimensions into the sum of two gauge invariant terms with one of them depending 
on the beta-function. The latter term provides conformal symmetry breaking contribution to the 
dilatation operator both in QCD and SYM theories with M = 1, 2. As follows from our analysis, 
the pairing of the eigenstates persists in QCD for terms depending on the beta function, but 
does not hold for the remaining "conformal" piece. 

Starting from two loops, the dilatation operator in SYM theories depends on the number of 
supercharges M . This dependence comes about through the contribution of 2(J\f — 1) real scalars 
and M gaugino fields propagating inside loops. We find that both contributions are irrelevant 
for enhancement or insuring of integrability since the M— dependence of the two-loop dilatation 
operator can be factored out (modulo an additive normalization factor) into a multiplicative 
c-number. This property makes the eigenspectrum of the two-loop dilatation operator alike in all 
gauge theories including the M = 4 SYM. Thus the main conclusion of our analysis is that the 
phenomenon of integrability is only tied to the planar limit and is sensitive neither to conformal 

4 The two charges commute for the operators with zero total light-cone momentum only, the so-called forward 
limit. 

5 One should mention however that, in general, integrability does not necessary imply degeneracy in the eigen- 
spectrum. The simplest example is provided by integrable open spin chains which determine the one-loop dilata- 
tion operator in multi-color QCD in the sector of mixed antiquark-gluon-quark operators ^|SJ- 
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symmetry nor supersymmetry. 

In our subsequent consideration we thoroughly elaborate on this subject. The presentation is 
organized as follows. In the following section we describe gauge theories under study — QCD and 
a generic Yang-Mills theory with fermions and scalars belonging to the adjoint representation 
of the gauge group, — and specify three-particle operators whose renormalization is scrutinized 
here. Then in Sect. El we work out two equivalent representations of the dilatation operator — 
as a mixing matrix on the space of local Wilson operators and as an integral operator in the 
representation of momentum fractions. The former representation is convenient for discussing 
symmetry properties whereas the latter one is the most appropriate for higher-loop calculations. 
To start with, we describe in Sect. 01 the main features of the approach on a simpler example 
of the one-loop dilatation operator. Then, we move on to the two-loop case and compute the 
irreducible components of the two-loop three-particle kernel term by term in Sect. El The reader, 
who is not interested in intermediate details but only wants to see the final result, is directed 
straightforwardly to Sect. 15.41 Having computed the two-loop dilatation operator, we turn to the 
study of its eigenspectrum in Sect. H3 At first, we address in Sect. 16.11 the sector of two-particle 
(twist-two) operators and diagonalize the dilatation operator at one loop making use of the 
conformal invariance. Then, we describe the constraints on the form of the dilatation operator 
imposed by conformal symmetry breaking in gauge theories. In Sect. 16. 1.51 we calculate the two- 
loop anomalous dimensions of twist-two operators in SYM theories with an arbitrary number of 
supersymmetries M and demonstrate that all of them are expressed in terms of the beta-function 
of the underlying gauge theory and the so-called "universal" anomalous dimension. Finally, 
we turn in Sect. 16.21 to the three-particle dilatation operator. We review the relation between 
integrability of the one-loop dilatation operator and pairing of its eigenstates and demonstrate 
that the double degeneracy observed in QCD and all SYM theories to one loop implies the 
existence of a conserved charge. We extend this consideration to two loops in Sect. 16.2.21 and 
explicitly construct the set of operators commuting with the dilatation operator and among 
themselves. We present the two-loop spectra of anomalous dimensions in Sect. 16.2.31 for QCD 
and all SYM theories and draw conclusions about two-loop integrability. We close this section 
with presenting two limiting cases in which we are able to reconstruct the all-loop dilatation 
operator and identify the corresponding integrable structures. Section [7| contains concluding 
remarks. Two appendices provide technical details of the calculational procedures used in our 
analysis. 

2 Case of study 

The goal of the present study is to reveal hidden integrability in four-dimensional Yang-Mills 
theory describing the coupling of gauge fields to fermions/scalars. Depending on the represen- 
tation in which the latter fields are defined one can distinguish two different types of the gauge 
theories: QCD and supersymmetric extensions of Yang-Mills theory. In the first case, the gauge 
fields are coupled to quarks in the fundamental representation of the SU(N C ) gauge group and 
the QCD Lagrangian is given by a well-known expression 

£qcd = -|trF^F^ + i$ lv p^ ■ (2.1) 

Here quarks are described by four- component Dirac fermions ip = (X a ,x a ) an d the gauge field 
strength F^ v = ^[D^D V ] is determined in terms of the covariant derivatives = 3^ — igA^t" 1 
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with generators t a in the fundamental representation of the SU(N C ) normalized conventionally as 
tr {t a t b ) = \5 ab . In the second case, the gauge fields are coupled to fermions and scalars belonging 
to the adjoint representation of the SU(N C ) group. The Lagrangian of the corresponding Yang- 
Mills theory takes the form 

£ adj = tr { - \F, U F^ + 2f\ &A af D"\f + \ {D^ AB ) (D^ AB ) 

- h Y \ aA $ AB , Af ] + h Y A^[0 AB , \%] - h 4 {0 CD , AB ] [0 AB , 4> CD \ } , (2.2) 

with the gauge fields, gauginos and scalars, X = {A^, X A , 4> AB }, all in the adjoint representation 
of the gauge group, X = X a t a and a,b = l,...,J\f. The covariant derivative in the adjoint 
representation is given by = <9 M — iglA^, }. The gauginos are described by the Weyl fermion 
X A which belongs to the fundamental representation of an internal SU(Af) symmetry group with 
its complex conjugate \ A = (\ A )*. The scalars are assembled into the antisymmetric tensor 
(p AB = —<p BA , with its complex conjugate (cp AB )* = 4> AB - The supersymmetric Yang-Mills 
theories with M = 1, 2 and 4 supercharges are obtained from the Lagrangian (j2.2j) by adjusting 
the values of the coupling constants and the number of gaugino and scalar species, rif and n s , 
respectively: 

• M = 1 SYM: there is a single specie of the gaugino, n/ = 1, and no scalars, n s = 0, so that 

hy = 0, /l 4 = 0. (2.3) 

• M = 2 SYM: there are two species of gaugino, rif = 2, and two real scalars, n s = 2, 

h Y = g, K = ±g\ ^ AB = V2e AB ^. (2.4) 

• M = 4 SYM: there are four species of gaugino rif = 4, and six real scalars, n s = 6, 

hy = V2g, h 4 = y 2 , <p AB = le ABCD ^ CD . (2.5) 



The above relations between the coupling constants are imposed by supersymmetry. It also 
relates to each other the number of gauginos and scalars 

Tl 

n / = l + y=AT. (2.6) 

One of our findings is that, as we will show in Sect. l5.2?2| integrability is not tied to supersymme- 
try. In other words, the phenomenon persists in the generic Yang- Mills theory ()2.2j) for arbitrary 
Af, to two loop order at least. 

The classical Lagrangians of the Yang-Mills theories (|2.1|) and (j2.2j) enjoy the conformal 
symmetry which will play an important role in our analysis. On the quantum level this symmetry 
is anomalous (except of Af = 4 SYM) due to running of the gauge coupling constant, 

dln 9 2 al 2\ a 9 2 a ( g 1 " 



The beta-function in QCD is given to two loop order by the well-known expression 

/WS 2 ) = ~£- 2 (yAT c - | n/ ) - ( (yiV c 2 - C F n f - 5 -N c n f ^ + <D(g«) , (2.8) 



where Cp = (N% — 1)/(2N C ) is the Casimir operator in the fundamental representation of the 
SU(N C ) group and n/ the number of the quark flavors. In supersymmetric Yang-Mills theory 
with M = 1,2,4 supercharges the two-loop beta-function is given by 

/W<? 2 ) = -^jr(4 -W- (0) 2 (2 - AO (4 - AT) + CV) . (2.9) 

For A/" = 4 the beta-function vanishes to all orders of perturbation theory and the gauge theory 
remains conformal on the quantum level. For M < 4 the beta-function is nonvanishing and, as 
a consequence, the conformal symmetry is broken. For Af = 2 the beta-function receives contri- 
bution to one-loop only while for M = 1 the beta-function coincides with the QCD expression 
(12. 8|) upon substitution rif — > N c and Cp — > N c . 

The central object of our study is the scale dependence of the gauge invariant nonlocal 
operators given by the product of elementary fields X{zn ll ) "living" on the light-ray defined by 
a light-like vector (so that n 2 ^ = 0). The explicit form of these operators depends on the 
representation of the SU(N C ) gauge group to which these fields belong. For elementary fields 
X = X a t a in the adjoint representation one can construct single-trace operators of the length L 

O adi (z 1 ,z 2 ,...,z L )=tx{X(z 1 n)X(z 2 n)---X(z L n)} , (2.10) 

where the gauge invariance is restored by inserting the Wilson lines in between each pair of fields, 

[zj,z j+ i] — Pexp ^ig J dzn^A^(zn)^j . (2.11) 

Later on we will adopt the light-light gauge (n ■ A) = A + = 0, so that the gauge links all shrink 
to the unit matrix. In what follows we do not display them for brevity. For fields X J (quarks) 
in the fundamental representation, we can introduce a baryon operator 

O fun (z x , z 2 , . . . , z Nc ) = e jlh ... jNc X^{z 1 n)X^(z 2 n) ■ ■■X^{z Nc n) . (2.12) 

It is built from exactly N c fields and is completely antisymmetric with respect to permutation of 
any pair of them. 

The nonlocal light-cone operators (|2.1()|) and (j2.12j) have the meaning of generating functions 
of local gauge invariant Wilson operators of the maximal Lorentz spin. The latter operators can 
be obtained from Taylor expansion of O a dj(^i, z 2 , ■ ■ ■ , Zl) and Of un (;?i, z 2 , . . . , zn c ) in powers of Zi. 
Each field entering ()2.10j) and ()2.12j) can be decomposed into different helicity components and, 
as a consequence, the nonlocal light-cone operators can be classified according to their helicity 
content. Among them one can distinguish light-cone operators with all helicities aligned. In 
what follows we shall refer to them as maximal helicity operators. 

It is well-known that local Wilson operators mix under renormalization and satisfy the evo- 
lution (Callan-Symanzik) equations. The same is true for the light-cone operators ()2.10j) and 
()2.12|) although the explicit form of the evolution kernels driving the scale dependence of these 
operators is different due to their nonlocal form. In particular, for nonlocal light-cone operators 
(J2.1(J|) built from the so-called "good" components of fundamental fields the evolution equation 
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takes the following form 6 j22l E3] 

U% + ^9 2 )-^j O adj (zi, z 2 , . . . , z L ) = -[H L ■ O adi \( Zl , 22, ... , z z ) + 0(l/iV c 2 ) , (2.13) 

where 0(1/N^) denotes terms corresponding to transitions involving splitting of the single trace 
into multi trace operators, i.e., O a aj{zi, £2, • • • , %l) — > ©adj(<2i, ^©adj^, . . . , z£). The evolution 
equation ()2.13|) follows from the Ward identity under dilatation transformations in the underlying 
gauge theory. 

The integral operator Hi entering (J2.13)) defines a representation of the dilatation operator D 
on the space spanned by nonlocal light-cone operators ([2. 10)1 . To the lowest order of perturbation 
theory, the explicit form of Hi in Yang- Mills theories with an arbitrary number of supercharges 
has been found in Ref. j23]. It turns out that to one-loop order, the evolution kernel has a number 
of remarkable properties in the multi-color limit: 

• The length L of the single trace operator is preserved and, as a consequence, Hi can be 
realized as a quantum mechanical Hamiltonian for a system of L particles with nearest- 
neighbor interaction, Hi = Hyi + if 23 + ■ ■ ■ + Hl\\ 

• Hi reveals a hidden integrability symmetry in the sector of maximal helicity operators 
both in QCD and in supersymmetric Yang-Mills theories: it possesses exactly L conserved 
charges and, therefore, is completely integrable; 

• Hi can be mapped into a Hamiltonian of the Heisenberg spin chain of length L and spin 
operators being generators of the SL(2\M) group with Af counting the number of super- 
charges. 

The main question that we would like to address in the present paper is whether the integrability 
phenomenon just described is preserved to two-loop order and beyond it. 

To this end, we shall evaluate two- loop corrections to the evolution kernel Hi and study its 
properties. Due to different partonic content in Yang-Mills theories with different number of 
supercharges, the number of maximal helicity light-cone operators varies with TV". Calculation of 
the two-loop evolution kernel for any of these operators is an extremely tedious task. Supersym- 
metry implies that the evolution kernels corresponding to different maximal helicity operators of 
the same length L are related to each other. Making use of this property, we shall restrict our 
analysis to maximal helicity operators built from the helicity +| component of the gaugino field 

(with a = 1, 2) 

At = K^ +/ ^=( A ^ . (2.14) 



2 U a(3» 7 \ 

It actually possesses a single nonvanishing component X A = X aA t a (with a = 1, ...,iV| — 1 
and a = 1, . . . ,J\f) which is charged with respect to internal SU(Af) isotopic group. Then, one 
decomposes the product of L fields over irreducible SU(J\f) components and selects the one with 
the maximal i?-charge 

®{z x ,z 2 ,...,z L ) = S tr {\ Al (zm) \ M (z 2 n) ...\ AL (z L n)} , (2.15) 

A 1 A2...A L 



6 This equation has been written under tacit assumption that the light-cone operators H2.1t)|l evolve au- 
tonomously in the multicolor limit. In general, Hl has a matrix form as the light-cone operators with different 
partonic content could mix with each other. 
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where the operation S stands for symmetrization over the 5 , f/(A/') indices. The counter-part 
of the operators 1)2.15)1 in QCD is the baryon operator (J2.12)) built from helicity ±| "good" 
component of the quark fields 

= K 1 -75)7-7+^ 






v y 



^+1 = 1 (! + 75)7-7+^ 



/0\ 




V 9/ 



(2.16) 



which have a single nonvanishing component transforming in the fundamental representation of 
SU(N C ). Later in the paper, we will encounter the operators of length L = 2 and L = 3: 

• For L = 2, the generating functions of maximal-helicity operators in SYM and QCD read, 
respectively, 

©adj(^i, z 2 ) = ti[X^ A ( Zl n)X B \z 2 n)] , (2.17) 
Ofim(zi, 22) = l^+l(z 1 n)'j + 'y ± ip + ^(z2n) = ^{z x n)q- >(z 2 n) , (2.18) 

where g J and are independent helicity ±| components of the quark Dirac field ^ + and 
7_i_ = (71 — 272) / y/2 is the antiholomorphic Dirac matrix. 

• For L = 3, analogous generating functions look like 

OadjOzi^Za) = S tr{\ A { Zl n)\ B {z 2 n)\ c {z 3 n)} , (2.19) 



ABC 



"fun 



zi,z 2 ,z 3 ) = e jk iq 3 (z 1 n)q k (z 2 n)q l (z 3 n) . 



(2.20) 



We remind that for fermions in the fundamental representation of the SU(N C ) group, the length 
of the operator ()2.12)) ought to be N c . Therefore, the operator Of un (zi, z 2 , 23) is well defined 
only for N c = 3. The operators Of UQ (zi, z 2 ) and Of un (zi, z 2 , z 3 ) have a direct phenomenological 
significance: their matrix elements determine the transversity distributions in the nucleon |2H] 
and the distribution amplitude of the delta-isobar [2E], respectively. 

We would like to stress that the choice of the maximal helicity operator is ambiguous. The 
reason why we prefer to work with the operators built from helicity +| fermions (gaugino, quarks) 
is that the number and complexity of contributing Feynman diagrams is gradually reduced and 
more importantly we can make use of available two-loop QCD calculations. In addition, as we 
will show in Sect. El the two- loop dilatation operator acts elastically on the maximal helicity 
operators built from the fermion fields ()2.14j) and ()2.1fij) . In other words, the operators ()2.15j) 
evolve autonomously and the corresponding evolution kernel ELj, can be realized as a quantum 
mechanical Hamiltonian. To two-loop order, H.l is given by the sum over two- and three-particle 
kernels, H^k+i and Hk,k+i,k+2, respectively. To identify the explicit form of these kernels it 
suffices to consider maximal helicity operators of length L = 3, Eqs. (|2.19jl and ()2.2())1 . Our 
strategy will be to evaluate the two-loop evolution kernels for the operators ()2.19j) and ()2.2()j) 
and, then, draw conclusions about their symmetry properties. 



3 Representations of the dilatation operator 

Before we proceed with actual calculations we would like to summarize general properties of the 
evolution kernels entering the renormalization group equation for the light-cone operators ()2.19j) 
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and (J2~20jl 

d d 

^ + P(g 2 )^ ) 0>(*i, z 2 , z 3 ) — — [H ■ 0](^, z 2 , z 3 ) , (3.1) 



where HI = Hl = 3. The evolution kernel EI admits the perturbative expansion 

e = a e (0) + A 2 e (1) + e>(A 3 ) , (3.2) 

with A = Cg 2 /87T 2 and the color factor C depending on the representation of the fermion fields. 
Its value will be specified in Sect. 14.21 

The central object of our study is the spectral problem for the dilatation operator 

W* q (Zi) = 7q(A) ■ (3.3) 

Here the eigenfunctions ^ q (zi) depend on the light-cone coordinates, 7 9 (A) define the eigenvalues 
of the dilatation operator and q denotes the complete set of quantum numbers to be specified 
later on. As we will show in Sect. 16. 2\ hidden integrability of the dilatation operator leads 
to a peculiar regularity in the spectrum of 7q(A). In gauge theories with unbroken conformal 
symmetry, (3(g 2 ) = 0, the solution to the evolution equation (|3.1j) looks like 

Q{ Zl , z 2 , z 3 ) = *2, 23) O,(0) , (3.4) 

Q 

with O q (0) being local Wilson operators. Substituting this relation into (j3.1|) one finds that 
the operators O q (0) have an autonomous scale dependence and their anomalous dimensions are 
given by 7 q (A). In gauge theories with fl(g 2 ) 7^ 0, the conformal anomaly induces the mixing 
between the operators O q (0) starting from two loops. It is driven by the additional term in 
the evolution equation (3{g 2 )d^ q {zi) / dg 2 proportional to the beta-function. It is straightforward 
to write down a solution to the evolution equation in that case but its explicit form is more 
complicated as compared to the one with (3(g 2 ) = 0. 

3.1 Mixing matrix 

We remind that the nonlocal light-cone operators (J2.19|) and (|2.2(J|) generate infinite towers of 
local Wilson operators. In particular, for the operator (j2.19|) one finds 

Q(z z z)~ V { ~ tZl)kl { ~ tZ2)k2 (-^ C (0) ( 35 ) 

k!,k 2 ,k 3 >0 l ' 2 ' 3 ' 

with the Wilson operators 

O klk2k3 (0)= S tr{(iD + ) k i\ A (0)(tD + ) k >\ B (0)(tD + ) k *\ c (0)} (3.6) 

ABC 

and D + = {n ■ D) = d + in the light-like axial gauge {n ■ A) = 0. Substituting (j3.5|) into the 
evolution equation (j3.1|) one recovers a conventional renormalization group equation for three- 
particle Wilson operators 

i^ + ^9 2 )^- 2 )o klk2k M = - E V El 3 <W0). (3.7) 
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Here the mixing matrix defines a representation of the dilatation operator ()3.2j) in the basis of 
polynomials \ki,k 2 ,k 3 ) = (—izi) kl (—iz 2 ) k2 (—iz 3 ) k3 /(ki\k 2 \k 3 \) 

M\h,l 2 ,l 3 )= ttl*i.*»'*>)' ( 3 - 8 ) 

To find the spectrum of the dilatation operator ()3.3|) one has to solve the spectral problem for 
the evolution kernel or equivalently diagonalize the mixing matrix 

E V S 3 Qifafc («) = 7,(A) c fclfc2fc3 (fl) • (3-9) 

2i,22,*3>0 

The coefficients Ck^k^q) define the expansion of the eigenstates ^ q {zj) over the basis ()3.8j) 

u^s) = e ( i z ; )ta ( i z f ( tf ■ (3.10) 

rrl ^2! fe! 

Let us take into account that the mixing can occur between the Wilson operators Ok 1 k 2 k 3 {fy with 
the same canonical dimension. This implies that the mixing matrix V^'j^. has nonvanishing 
entries only for k\ + k 2 + k 3 = li + l 2 + h- in terms of the evolution kernel, Eq. (j3.8j) , this property 
states that the evolution kernel preserves the overall power of the polynomials \ki, k 2 , k 3 ), or 
equivalently 

[M, z,d zi + z 2 d Z2 + z 3 d Z3 ] = . (3.11) 

In addition, due to Poincare covariance, the anomalous dimension of the operators Ok^k^ix) 
does not depend on the coordinate x. In terms of the evolution kernel, this implies that EI has 
to be translation invariant 

[H, d Zl + d Z2 + d Z3 ] = . (3.12) 

Since the dilatations and translations do not commute with each other, the eigenstates can not 
diagonalize the two charges simultaneously unless the total light-cone momentum equals zero, 
(d zl +d Z2 +d Z3 )y q ( Zl ) = 0. 

Putting z\ = z 2 = z 3 = in (J3.5)) one obtains a local operator with no derivatives 

0(0,0,0) = O 000 (0)= l -f abc S X aA (0)X bB (0)X cC (0). (3.13) 

4 ABC 

This operator evolves autonomously and according to (|3.7j) and ()3.8|) its anomalous dimension 
defines the eigenvalue of the dilatation operator 

e|o,o,o) = 3r(A)|o,o,o) , (3.14) 

where the additional factor 3 was inserted for the later convenience. As we will show in Sect. 16.273*1 
the anomalous dimension 

r(A) = Ar (0) + A 2 r (1) + c(A 3 ) (3.15) 

sets up the lower /upper bound in the spectrum of three-particle dilatation operators (|3.2[) . 
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3.2 Momentum representation 

The evolution operator H acts along the light-cone direction in Minkowski space-time. Yet an- 
other representation of EI can be obtained by going over from the configuration to the reciprocal, 
momentum space. To this end, one performs Fourier transformation of the light-cone operators 

/oo poo poo 

d Ul du 2 du3^ UlZl+iUzZ2+iU3Z3 0{ Ul ,u 2 ,u 3 ) (3.16) 
-oo J — oo J — oo 

with Ui having the meaning of light-cone momenta of particles. Substituting ()3.5j) into this 
relation one finds that the local Wilson operators are given by the moments of 0(ui, u 2 , u 3 ) 

O klk2k M = j duxd^dusU^u^u^Oim,^,^). (3.17) 

In the momentum representation, the Callan-Symanzik equation is known in QCD as the Brod- 
sky-Lepage equation 



d d \ ~ f ~ 

^dfi + ^^di?) <D> ^ 1 ' U2) = J ^ 3V ( Ml ' U2 ' u ^ Vl ' U2 ' ^3)0(^1, v 2 , v 3 ) , (3.18) 

where the integration measure is 

/ 3 3 \ 

[dv} 3 = dv\dv 2 dv 3 5 I ^Jfj — ^^Wj J • (3.19) 
\i=i 3=1 J 

The evolution kernel in the momentum representation, V, is obtained from EI through the Fourier 
transform 

H. e - iS ^ fc2fc = - J [du] 3 Y(u\v) e-^ UhZk , (3.20) 

with [du} 3 determined by Eq. ()3.19|) with U{ and V{ interchanged. According to (|3.12|) . the 
evolution kernel preserves the total momentum v 1 + v 2 + v 3 = u\ + u 2 + u 3 . This property is 
automatically preserved in (|3.18|) . In addition, it follows from (|3.11|) that the evolution kernel is 
a homogeneous function of the momentum fractions 

V(Au 1 , Xu 2 , Xu 3 \Xv l , Xv 2 , Xv 3 ) = A~ 2 V(ui, u 2 , u 3 \v 1: v 2 , v 3 ) . (3-21) 

Obviously, the perturbative structure of the evolution kernel V is identical to the one introduced 
for H, Eq. (13~2|) . 

Combining together (J3.2(Jj) and (|3.8j) one obtains the following representation for the mixing 
matrix 

f[du] 3 ¥(u\v)u k M 2 u k 3 3 = - £ <[t\A^i (3-22) 

J h,h*h>° 

that is, positive integer moments of the evolution kernel with respect to u— variables are polyno- 
mials in v— variables of the same degree k\ + k 2 + k 3 = li + l 2 + l 3 . Below this property will play 
an important role in our analysis and we shall refer to it as the polynomiality condition. 
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By construction, the eigenvalues of the evolution kernel Y(u\v) defined in 1)3.20)1 coincide with 
7q(A), Eq. (J3.3)) . To see this, one considers the spectral problem for the mixing matrix 

™ klk2hi (Q) V^Sa = 7 fl (A) ™ hhh (q) . (3.23) 

Comparing this relation with ()3.9|) one observes that the coefficients w lll2la (q) and Ck^k^q) are, 
respectively, the left and right eigenstates of the same mixing matrix corresponding to the same 
eigenvalue 7 q (A). Introducing the polynomial 

P q (ui) = Yl ^^{q^u^ul 3 (3.24) 

k\,k2,kz>0 

one finds from ()3.22|) that it diagonalizes the evolution kernel 

J [du} 3 Y(u\v) P q { Ui ) = - 7 ,(A) P q {vi) . 

The left and right eigenstates of the mixing matrix are orthogonal to each other 
in terms of the polynomials (j3.1()j) and ()3.24|) this condition reads 

Pq{d zl ,d Z2 ,d z:i )^ q >{z 1 , z 2 , z 3 )\ z _ Q ~ 5 qql . (3.26) 

Combining this relation together with ()3.4|) one concludes that the P— polynomials determine 
the explicit form of local Wilson operators entering (|3.4J) 

O q (0) = P q (d Zl ,d Z2 ,d Z3 )O(z u z 2 ,z 3 )\ Zi=0 . (3.27) 

Moreover, it follows from ()3.21|) that the eigenstates of the evolution kernel V(u\v) have to be 
homogeneous polynomials 

P q {\ui, Xu 2 , Xu 3 ) = X N P q (u 1 , u 2 , u 3 ) . (3.28) 

Here, integer N > counts the total number of light-cone derivatives in the expression for the 
Wilson operator (j3.27J) . It also defines the Lorentz spin of O q (0). The mixing between the 
operators with different N is protected by Lorentz symmetry and, therefore, N can be identified 
as one of the quantum numbers q. 

In this section we described the dilatation operator in two different representations — as a 
mixing matrix for local Wilson operators and as an integral operator in the representation of 
momentum fractions. Obviously, the spectrum of the dilatation operator does not depend on the 
particular representation and its choice is a question of convenience. As we will demonstrate in 
Sect. 03 the momentum representation has a number of advantages as far as the two-loop calcula- 
tion of the evolution kernels is concerned. Firstly, by virtue of Lorentz covariance, the eigenvalues 
of the evolution kernel W(u\v) do not depend on the total light-cone momenta Ylj u j = Tlij v j- 
This allows one to put its value equal to zero. From the point of view of nonlocal operators, this 
corresponds to assuming translation invariance of Q(zi, z 2 , . . . , zl) along the light-cone, or, equiv- 
alently, to considering forward matrix elements of the nonlocal light-cone operators in (J3.16|) . In 
this way, one automatically eliminates the contribution of local operators with total derivatives 
in ()3.5j) . which in turn simplifies significantly the form of the mixing matrix (|3.7j) . Secondly, 
there exist well developed techniques for calculating the renormalization group kernels in the 
momentum representation. We shall profit from higher-order results available in the literature. 



(3.25) 
. Being written 
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4 Dilatation operator: 1-loop 



Let us now turn to the computation of the evolution kernel H for the maximal helicity operators 
(j2.19j) and ()2.20|) both in QCD and supersymmetric Yang-Mills theory. To start with, we review 
in this section the one-loop calculation of H. It will serve as an illustration of the general 
formalism which will then be used for more involved two-loop analysis. 

4.1 Regularization scheme 

Performing the calculation we shall adopt the light-like axial gauge (n-A(x)) = A + (x) = 0. Since 
the evolution operator H is a gauge-invariant quantity, any choice of gauge fixing is possible. In 
covariant gauges, the nonlocal light-cone operators ()2.1()j) involve additional gauge links (|2.11j) 
between the elementary constituent fields. This leads to proliferation of the number of Feynman 
diagrams contributing to EI. The light-like axial gauge A + {x) = is advantageous as the number 
of relevant graphs is much smaller, especially at higher loops. Yet another advantage of this 
choice is that we can employ certain results for one- and two-loop evolution kernels in QCD 
available in the literature. 

The Yang-Mills theory in the light-like axial gauge has a number of subtleties j2Z|. The gauge 
field propagator (0|T A a {x) A b u (y)\0) = (—i)S D IJiV (x—y) is given in the momentum representation 
by 

^ )= PTl' d^(k) = g, u - ^ ^ KU " (4.1) 

with k + = {k-n). It has a spurious singularity at k + = which indicates that the gauge ambiguity 
is not fixed completely for the choice A + {x) = 0. To give the meaning to the gluon propagator 
(|4.1j) . one has to specify the prescription for integrating around the pole k + = inside momentum 
Feynman integrals. The only available residual gauge fixing which is consistent with causality 
properties of Feynman integrals is the one due to Mandelstam and Leibbrandt [211 [2H] 



[fc+] ML k+k_ + iQ- 

It has been checked for various gauge invariant quantities that to two-loop accuracy this pre- 
scription leads to the same results as the Feynman gauge. However, the calculations within the 
Mandelstam-Leibbrandt prescription are rather involved and there exist another popular choice, 
the so-called principal value prescription, defined as [^01 

1 1/1 1 

[fc+] pv 2^+ + ^ k + - i5 

with 5 —>■ 0. This prescription has proved to be viable and robust in higher order QCD calcula- 
tions of the evolution kernels and we will accept it for two-loop calculations in our subsequent 
analysis. 

To isolate divergences of Feynman integrals we will employ in parallel two different regulariza- 
tion schemes — dimensional reduction (DRED) and dimensional regularization (DREG). In both 
schemes one regularizes Feynman integrals by setting the space-time dimension to D = 4 — 2e 
and taking the limit e — > afterwards. In the DRED scheme one keeps the number of all tensor 
fields components (gauge bosons and fermions) to be the same as in four dimensions (HI] while 
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in the DREG scheme this number is e— dependent. The two schemes are equally suitable for 
two-loop calculations but the usage of the DRED is mandatory as long as one wants to preserve 
the supersymmetry of Yang-Mills theory. 

To determine the evolution kernel H we shall examine perturbative corrections to the matrix 
elements of nonlocal light-cone operators ()2.19|) and ([2.20)1 . For D^4, the ultraviolet divergences 
manifest themselves as poles in e and we shall renormalize "bare" nonlocal light-cone operators 
using the modified minimal subtraction procedure. Depending on the choice of the DREG or 
DRED regularization schemes, it will be denoted as MS and DR, respectively. The renormalized 
operator is defined as 

0*0*1, z 2 , z 3 ) = Z- 0< bar< %, z 2 , z 3 ) , Z(l/e, g 2 ) = l + J2 r 1 ' ^ 

n=l 8 

with the subtraction "constants" Z^(g 2 ) chosen in such a manner that the Green functions 
(0|O R (2;i, z 2 , Z3) .. . |0) remain finite for e — > 0. 7 Substituting ()4.4|) into the evolution equation 
f[2.13|) one finds that the evolution operator is related to the residue of the simple 1/e pole 

(ilnZfl/e, q 2 ) d m, 9x 

M = ,\ 1 = -: — Z [1] (g 2 ). (4.5) 

din// ding v ' K ' 

Here in the second relation we used the fact that the gauge coupling constant acquires a nontrivial 
dimension for e 7^ and, as a consequence, the beta-function in the D-dimensional Yang-Mills 
theory possesses an additional e— dependent term 

^ = p £ (g 2 ) = -2e + P(g 2 ). (4.6) 

As we will explain in detail in section lfi.l.H[ the presence of this term has far reaching consequences 
for the properties of the dilatation operator. In particular, it implies that in Yang-Mills theory 
with f3(g 2 ) = the conformal symmetry is broken for e 7^ 0. 

We would like to stress that the renormalization constants in the MS— and DR— renormali- 
zation schemes are different but they are related to each other though a finite renormalization. 
As a consequence, the evolution operator depends on the renormalization scheme and one has 
to distinguish the operators H^s and M^r. As we will show in Sect. 15.4.11 these operators differ 
starting from two-loop order only and the difference amounts to a c-number. As a result, if 
the dilatation operator possesses a hidden symmetry it will hold independently of the chosen 
renormalization procedure. 



4.2 One-loop 

Examining Feynman diagrams contributing to 0(zi, z 2 , z 3 ) it is easy to see that to the k th order 
in A the interaction could occur between (k + 1) particles at most. As a consequence, to one- 
loop order rf ) receives contribution from pair-wise interactions as well as their self-energies. To 
determine the one-loop evolution kernel, one has to calculate one-loop corrections to the nonlocal 
light-cone operator 0(z\, z 2 , z 3 ), renormalize it according to (|4.4j) and, finally, evaluate H (0) using 

7 It worth mentioning that the renormalization constants entering H4.4|l are certain integral operators acting on 
the light-cone coordinates of < - barc ^ ) (zi, Z2, Z3). 
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(a) (b) (c) 

Figure 1: One- loop diagrams in the light-cone gauge contributing to the pair- wise kernel (a) 
and the self-energy due to gluons (6) and scalars (c). 

(14. 5 j) . We sketch below the intermediate steps while the details of the calculation can be found 
in Refs. 

Since to one-loop order the interaction can only occur between two fermions, the dilatation 
operator has a pair-wise structure 

H (0) = + + + 3rW . (4.7) 

Here the two-particle kernel M.^ acts locally on the j th and k th fields while 1^°) receives contri- 
bution from the self-energy corrections to the fermion. They are determined by the Feynman 
diagrams shown in Fig. ^ Notice that in the SYM theory, the self-energy diagrams involve both 
gauge field and scalar exchanges while the diagram for contains solely the gluon exchange. 
The reason for this is the following. One can straightforwardly deduce from the Lagrangian (|2.2j) 
that the scalar exchange corresponds to the transition \* A \ B } — > Xc^d (symmetrized over a and 
b) involving the scalar propag ator (0 AC BD ). As follows from (|2~5jl . this propagator vanishes 
for Af = 1 and Af = 2, while for Af = 4 it is proportional to ~ £ abcd an j vanishes upon 
symmetrization over the isotopic indices a and b. 

The contribution of the diagram in Fig. da) to the one-loop dilatation operator in QCD and 
SYM only differs by an overall color factor. For gauginos in the adjoint representation, the color 
factor equals T^Tg with (Tj)% c = %f ahc being the SU(N C ) generators in the adjoint representation 
acting on the color indices of the j th particle. Since the three-fermion operators have zero total 
color charge, one has Y^j=i 2 3 ^7 = an d, therefore, 

T«n = l -{T?y - \{T«y - \{ny = -\n c , m 

with Ti = —{Tj + Tk) a and (Tj 1 ) 2 = N c being the Casimir operator in the adjoint representation 
of the SU(N C ). For quarks in the fundamental representation of the SU(N C ), the color factor 
equals 

£ i 1 ...i' j ...i' k ...i Nc {t )i' k ik = £ h---ij ■■■i k ---iN c ( J ' (^-9) 

where one takes into account the identity (t a )jk(t a ) nm = 8j m 8kn/2—8jk5nm/ (2iV c ). We remind that 
the baryon operator (|2.2()jl is only well-defined for N c = 3 in which case —(N c + 1)/(2N C ) = —2/3. 
As a result, the contribution of the diagram in Fig. ^a) to the one-loop evolution kernel for the 
light-cone operators (j2.19J) and (j2.20j) is accompanied by the prefactors 

A fun = J^-, A ^j = |^ (4-10) 
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correspondingly. In what follows we shall use Af un and A a dj as parameters of the perturbative 
expansion of the dilatation operator (|3.2j) . In this case, the two-particle kernel HI^? takes the 
same form for the light-cone operators in the adjoint, Eq. (j2.19j) . and fundamental, Eq. ([2.20)1 . 
representations 



Zi,z 2 ,z 3 ) = I — a 2j 1 



da 



o 



a 



2®( Zl ,z 2 ,z 3 ) 



[azi + az 2 , z 2 , z 3 ) 



0>(zi, azi + az 2 , z 3 ) 

(4.11) 

where a = 1 — a and j — 1 is the conformal spin of the quark/gaugino. For the later use, 
we shall display the dependence of the evolution kernels on the conformal spin j. According to 



Eq. ([4. II)) . the operator has a transparent physical interpretation: it displaces n ttl and k 
particles along the light-cone in the direction of each other. 

The gauge field contribution to the self-energy shown in Fig. ^b) is proportional to the 
quadratic Casimir operator: Ca = N c for the gaugino in the adjoint representation of the SU(N C ) 
and Cp = (N% — 1)/ (2iV c ) =4/3 for the quark in the fundamental representation of the SU (N c = 
3) group. In addition, in the SYM theory one has to include the contribution of scalars to gaugino 
self-energy (see in Fig. IHc)) which is proportional to ~ n s = 2(Af — 1). Combining together 
the two contributions one obtains the following expressions for the normalization constant T^ ) 
in (HZZD 

-(o) _ 1 r (o) I 



th 



QCD 



SYM 



-J\f. 



(4.12) 



Notice that T^ym coincides with the QCD expression for J\f = 1. Combining together ()4. 7j) . (|4.11|) 
and (j4.12j) . one obtains the one-loop dilatation operator for the three-particle quark and gaugino 
operators in QCD and SYM theories, respectively, in the coordinate representation. As follows 
from ()4.11j) . the two-particle kernel annihilates the local operator O(0, 0, 0) and, therefore, 3Fsy M 
determines the anomalous dimension of the local gaugino operator ()3.13j) . in agreement with 

dsm. 

Let us transform the obtained expression for IH into the momentum representation. Ap- 
plying ()3.20j) and performing the Fourier transformation one finds that Y(u\v) has the same 
structure as (|4.7j) 

V<®(u\v)= [Y ( -°\u 1 ,u 2 \v 1 ,v 2 )] + (4.13) 
+ [¥^(u 2 , u 3 \v 2 ,v 3 )] + [Y^(u 3 , Ul \v 3 , Vl )] - 3T^5( Ul - v x )b{u 2 - v 2 ) , 



with the two-particle kernel given by a well-known expression 

2j_1 6(wi,t>i) , fu 2 \ 2j ~ l Q(u 2 ,v 2 



v(°)( Ml; 



U 2 \Vi,V 2 



M2Y 

v 2 ) 



Vx-Ux \V 2 J v 2 -u 2 
Here the notation was introduced for a generalized "step" function 



5(u\ + u 2 — v 1 — v 2 ) . 



Q(u n , v n ) = 6(u n )6{v n - u n ) - 0(-u n )6{u n - v n ) . 
The symbol [•••]+ stands for the plus-distribution 



[ip(u, v)] + = ip(u, v) — 5(u — v ) / du'(p(u', v) 



(4.14) 



(4.15) 
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with <p(u, v) being a test function. This distribution regularizes the end-point singularity of (j4.14j) 
for Vi — > Ui in such a way that the evolution kernel (J4.13)) has finite moments (|3.22j) . Notice that 
due to the total momentum conservation, = ^f n , Eq. ([3.19)1 . the two-particle evolution 

kernel ()4.14|) preserves the momentum of the third particle, u 3 = t> 3 . By the same token, the 
constant term in the evolution kernel (|4.13jl preserves the momenta of all particles, u n = v n with 
n = 1,2,3. 

The variables Vi and Ui have the meaning of the light-cone momenta of particles. Calculating 
the moments of the evolution kernel ()3.22j1 . one can assign to Vi arbitrary real values. For given 
Vi, the integration region over the u— variables is not arbitrary and it is determined by the step 
functions (|4.14|) . It is instructive to consider a special case, the so-called Brodsky-Lepage (BL) 
limit HI], 

< Vi, v 2 , v 3 < 1 , V i + v 2 + v 3 = 1 . 



It is easy to see that the step function reduces in this limit to 6(wj, Vi 
a consequence, 



BL 



6( Ui )0(vi 



BL 



2J-1 



{Ui)0(vi - Ui) 



+ 



Ui 

V\) r, - //| 

u 2 x 



(4.16) 
Ui) and, as 

(4.17) 



»2 



9(u 2 )9(v 2 - u 2 ) 
v 2 - u 2 



5{ui + u 2 - vi - v 2 ) . 



Substituting this expression into ((4.13)1 . one verifies that the possible values of the u— variables 
are restricted to the same simplex as in ()4.1(ij) 



< Mi, U 2 , U 3 < 1 



U\ + u 2 + u 3 = 1 . 



(4.18) 



This property has a simple physical interpretation. It implies that for all three particles moving 
forward in time along the light-cone direction n^, the interaction described by the evolution 
kernel ()4.13|) can turn none of the particles to propagate backward in time. Remarkably enough, 
this property is general enough to hold to all order of perturbation theory p!4*l 133] . Below in 
Sect. El we will demonstrate its validity to two loops. In the following exposition we will stick 
to the BL-representation since, on the one hand, it simplifies the step-function structure of the 
momentum space dilatation operator and, on the other hand, suffices for evaluation of anomalous 
dimension mixing matrix. A simple procedure [3^] for reconstruction of the complete evolution 
operator without the restriction to the simplex ()4.16|) is given in Appendix El 



5 Dilatation operator: 2-loops 

In the previous section we have calculated the one-loop evolution kernel for the three-particle 
gaugino/quark operators (|2.19|) and (|2.20|) . We demonstrated that it is given by the same, uni- 
versal expression ()4.11|) both in QCD and in SYM theories with arbitrary number of supercharges 
N '. The fact that the fermions are defined in different representations of the gauge group mani- 
fests itself in the different color factors accompanying the coupling constant, Eq. ()4.1Ujl . Let us 
now extend our analysis to two loops. As we will see, the above mentioned universality gets lost 
beyond leading order and one has to consider separately the gaugino and quark operators. In the 
former case, we will start with the M = 1 SYM theory and will gradually increase the diversity 
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of the particle's content by going over to M = 2 and M = 4 theories. In this way, we will be able 
to clearly identify whether emerging integrability phenomena are sensitive to the representation 
in which the fermions fields are defined as well as to the presence of scalars and, on top of it, 
sup er symmet ry. 

Our analysis will be entirely formulated in the momentum space, i.e., we will calculate the 
two- loop kernel V of the evolution equation (|3.18|l . It is straightforward to translate the obtained 
expressions from the momentum to the configuration space by applying ()3.2()j) . As we will argue, 
to two-loop order, the three-fermion operators ()2.19j) and ()2.2())1 evolve autonomously and the 
perturbative expansion of the evolution kernel reads in the momentum representation 

V(«|«) = \¥ {0) {u\v) + A 2 V (1) (it|t>) + C(A 3 ) , (5.1) 

where depending on the representation of the fermion fields, the coupling constant A is given 
by the corresponding expression in (|4.1()jl . The one- loop result was given earlier in Eq. ()4.13j) . 
In addition, to simplify the calculation we shall apply the Brodsky-Lepage limit and assume 
that the v— and u— variables take nonnegative values and belong to the simplices (14.16)1 and 
(|4.18|) . respectively. This does not lead to a loss of generality since the resulting expression 
for V(ii|i>) can be uniquely continued to arbitrary values of the scaling variables. As before, 
the calculation will be performed in the light-like axial gauge A + (x) = within the DR— and 
MS— renormalization schemes. 

The two-loop Feynman diagrams contributing to in the light-cone gauge can be separated 
into three different sets sorted according to the number of fermions involved in the interaction, 
i.e., one-, two- and three-particle irreducible contributions. Their contribution to the evolution 
kernel can be written as 

Y {1 \u\v) = Vg, + V« + V« + VS 1 ^ + Vg + V£ } + VS 1} + V?> + V« , (5.2) 

where the additional terms have been added to ensure the invariance of the evolution kernel 
under cyclic permutations of quarks. Denoting the kernel of the operators and as 
W^(ui,u 2 ,u 3 \v 1 ,v 2 ,v 3 ) and VW(«i, u 2 \vi, v 2 ), respectively, one has 

= ^ X \ui,Uj,u k \vi,Vj,v k ) , Vfjj =¥ {1) (u j ,u k \v j ,v k ) . (5.3) 

The three-particle irreducible kernel V^3 receives contributions from the diagrams shown in 
Fig. |21 The set of two-loop diagrams contributing to the two-particle kernel is displayed 
in Figs. 0] and El Among these diagrams one can distinguish those without scalars (see Fig. 0J). 
They are exactly the same as involved in renormalization of twist-two maximal-helicity operators 
in QCD IHHj • However, as we will see momentarily, the important difference with the twist- 
two operators is that for three-fermion operators (|2.19|) and ()2.20|) each pair of particles has a 
nonvanishing color charge and this affects the color factors accompanying individual diagrams. 
Finally, the self-energy corrections to gaugino/quark fields induce c-number corrections V, ; , see 
Fig. which have the structure 

=& ) 5(u 1 - v 1 )5(u 2 -v 2 ), (5.4) 

with £^ being the residue of the renormalization constant of i th fermion line to two loop order. 
Let us consider three different contributions to Y^'(u\v) one after another. 
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12 3 (a) (b) 



Figure 2: The three-particle contributions (with mirror symmetrical graphs tacitly implied) to 
the three-particle dilatation operator. 



5.1 Three-particle contributions 

The three-particle irreducible kernel receives contribution from Feynman diagrams which 
either represent an iteration of a single-particle exchange, or contain a triple gluon vertex (see 
Figs. El (a) and (b), respectively). In QCD and M = 1 SYM theory, such kind of diagrams can 
only involve gluon exchanges while in the M = 2 and M = 4 theories scalars may, in principle, 
propagate instead of gluons in internal lines. As it was explained in Sect. 14.21 by virtue of the 
maximal R— charge of the three-particle operator (j2.19|) . the scalar fields do not contribute to the 
two-particle kernel to one loop order. The very same argument applies to the three-particle 
irreducible kernel V^- In other words, the diagrams with scalar exchanges vanish identically 
and the evolution kernel takes the same form in SYM theories with M = 1,2,4. 

Another interesting feature of the kernel V^3 * s that the diagram in Fig. El (b) involving the 
three-gluon vertex vanishes identically. This property holds true independently on the represen- 
tation in which the fermion fields are defined and, therefore, it is valid both in QCD and SYM 
theories. Namely, for the quarks in the fundamental representation of the SU(3) the color factor 
C corresponding to the diagram in Fig. El (b) equals 

The two sides of this relation have an opposite parity under permutations j\ ^ j'2 leading to 
C = 0. In the similar manner, for the gaugino in the adjoint representation of the SU(N C ) the 
color factor equals 

jabc T a T b T c = _jabc T a T b^ T c + j>c^ = q ? (5 g) 

where (Tj')bc = if abc is the color charge of j th particle and f abc = —f acb are the SU(N C ) structure 
constants. Here in the first relation we used the color neutrality of the single-trace three-particle 
operator, i.e., Ej=i, 2 ,3 T / = °> and applied the identity f abc T^T b = -\N c T c y In arriving at (pTKJl . 
it was crucial that the color-singlet operator (J2.19)) is built from L = 3 particles. For L > 3, the 
contribution of the diagram in Fig. |21 (b) to the multi-particle single-trace operator is different 
from zero but it is suppressed in the multi-color limit as ~ 1/N% since the diagram is nonplanar. 

We conclude that the three-particle irreducible kernel V^3 i s determined by the Feynman 
diagram shown in Fig. El (a) and its mirror image. The two diagrams describe iteration of one- 
gluon exchange and the corresponding color factor is merely the square of the one-loop color 
factor. For gluinos in the adjoint representation it yields T^T^T^T^ = (— 2 ) > where we have 
used the identity ()4.8|) twice. For fundamental quarks, one uses in the same fashion Eq. ()4.9|) . 
In both cases the color factors can be absorbed into the redefinition of the coupling constants 
(J4.10|) and, as a result, has the same form in QCD and in all SYM theories. 
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The evolution kernel ^\2 3 {u\v) depends on two sets of momenta. As was already mentioned, 
we restrict our analysis to the values of variables belonging to the simplices ([4.16]) and ([4.18]) . For 
given values of the v— variables and ui+u 2 +u 3 = 1, the three-particle irreducible kernel Vi2 3 (i*|?j) 
has support inside the triangle < U\, u 3 < 1 and u\ + u 3 < 1 shown in Fig. EJ Three sides of 
the triangle corresponds to u% = 0, U2 = and u 3 = 0. Due to a larger variety of possibilities 
of momentum redistributions between three interacting fields, Vi2 3 (u|v) has a more complex 
structure as compared with the two-particle kernel (j4.17|) . A simple analysis demonstrates that 
there exist six distinct regions in the phase space displayed in Fig. El the diagram in Fig. El (a) 
contains two internal gluon lines carrying the light-cone momenta v\ — u\ and u 3 — v 3 and the 
virtual fermion line with the momentum V\-\-v 2 —Ui. For the mirror symmetric diagram the latter 
momentum reads v 3 + v 2 — u 3 . In distinction with the external u— variables, the momenta flowing 
through internal lines are not necessarily positive in the Brodsky-Lepage limit. For v\ > u\ we 
can readily identify two regions: region 3 for v 3 > u 3 and combined regions 1 and 2 for v 3 < u 3 . 
In these kinematical regions the internal fermionic line does not bring in any new restrictions 
since the momentum flow is always positive V\ + v 2 — 1*1 > 0. Considering now the situation 
of the negative momentum flow in the first gluon line, i.e., v% < u\, one immediately finds that 
the virtual quark/gluino line can have either positive (v\ + v 2 > «i) or negative (t>i + v 2 < u 2 ) 
momentum flow, which is laid over the restriction from the second gluon line V3 ^ U3. From 
these we find three regions in the phase space labelled as 4, 5 and 6. Notice that the boundary 
Vi+v 2 = U\ is reflected from the boundary u 2 = of the support region into the boundary w 3 = t> 3 . 
The consideration of the mirror symmetric diagram to Fig. |21 (a) allows one to further separate 
the regions 1 and 2 from the positivity/negativity of the momentum flow of the internal fermion 
line. The above analysis covers all regions in the support region of the two-loop three-particle 
irreducible kernel. 

Summarizing, one decomposes the three-particle irreducible kernel as follows 

V^ 3 (u\v) = 9(u 3 - vi)¥x(u\v) + 9(v x - u 3 )9{v 1 - Ul )9(u 3 - ^ 3 )F 2 (^|^) (5.7) 
+ 9( Vl - ui)9(v 3 - u 3 )¥ 3 (u\v) + 9(u x - v x )9(u 3 - v 3 )¥ A (u\v) 
+ 6(v 3 - ui)F 5 (u|w) + 9(u! - v 3 )9(u! - v x )9{v 3 - u 3 )¥ 6 (u\v) , 

where Uj = 1 — Uj, Vj = 1 — Vj and the function ¥i(u\v) denotes the contribution of the i th region 
on the phase diagram in Fig. El It is tacitly assumed that the u— and v— momenta belong to the 
simplices (|4.18|) and ([4.16]) . respectively. The sum of the diagrams shown in Fig. |2]is symmetric 
under the interchange of the 1 st and 3 rd lines. Therefore the kernel (15.7]) has to be invariant under 
the interchange of the corresponding momenta {u\, u 3 \v\, v 3 } <-» {u 3 ,ui\v 3 ,vi}, which yields the 
following relations 

¥ 3 (u 1 ,u 2 ,u 3 \v 1 ,v 2 ,v 3 ) = ¥ 3 (u 3 ,u 2 ,u 1 \v 3 ,v 2 ,v 1 ) , (5.8) 
¥ 4 (u 1 ,u 2 ,u 3 \vi,v 2 ,v 3 ) = ¥ A (u 3 ,u 2 ,ui\v 3 ,v 2 ,v 1 ) , 
¥ 5 (ux,u 2 ,u 3 \v 1 ,v 2 ,v 3 ) = ¥ 2 (u 3 ,u 2 ,ux\v 3 ,v 2 ,v 1 ), (5.9) 
¥ e (u 1 ,u 2 ,u 3 \v 1 ,v 2 ,v 3 ) = ¥ 1 (u 3 ,u 2 ,u 1 \v 3 ,v 2 ,v 1 ) , 

so that the three-particle kernel ([5.7]) only involves four nontrivial functions. Moreover, as we 
will see in a moment, one of these functions vanishes. 

To determine the functions Fj one has to evaluate the three-particle connected diagrams 
displayed in Fig. El using particular renormalization scheme to treat ultraviolet divergences and 
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Figure 3: Phase- space regions of three-particle irreducible kernel. 

extract V^3(' u l t ') from a single pole ~ 1/e contribution using (j4.5J) . The calculation can be 
straightforwardly performed either in the covariant or light-cone formalism. We found that the 
contribution of diagrams in Fig. El to Y^i^-lv) is the same in the DR— and MS— schemes. It 
leads to the following expressions for the F— functions 

(ui - vi)(u 3 - v 3 ) 
{vi - u 3 ){u x - v 3 ) 



F i = 2J v vf( v ± - M 3, v 2 ) In — r - T7 z — (5.10) 



lr f \ ff - M Ml Ml -^3 .Iff - w,- \ i "l«l-"3 

+ o/K^li/K"! - ^3jhl^— r + of[U2,Vl -u 3 )}{vi - u 3 , v 2 ) In — - 

M2M 3 (« 3 -M 3 ) ^3(^1-^1 



TO \ ff— M ~ ^3) / K11 ^ 

F 2 = 5/(«i, vi)f{vi ~ u 3 , v 2 ) In — - (5.11) 

i f( srf - \ i Mi(^i-^ 3 ) 2 , 1,/ - \ f f — m ^(wi-va) 

+ 2j{ui,vi)f{u 2 ,ui - u 3 )m — — + 2/(^2, Mi - ^3)/(mi - ^3,^2) In = , 

U2M 3 (M3-M 3 ) M 2 W 3 

F 3 = 0, (5.12) 

w iff - \ ff- \ 1 ^1^3(^1 -%)(mi -^3) / K1 o\ 

F 4 = 2/(1*2, Mi - v 3 )f(ui - v 3 , v 2 ) In 2— (5.13) 

W|UlU 3 

Iff Sff - m Ml^2 i,, v,,_ M 3 (Mi-fi) 

+ 2/(mi, Wl)/(M 2 ,Ml - w 3 )ln — + 2/(Mi,fi)/(^i - M 3 ,f 2 )ln- 



^3lMi-fiJ Ui[Vi — U3) 

Iff - Sff m M3^2 ,1,/- N,/ M Wl(W3-«s) 

+ 2 j{u 2l V x -M 3 )/(M 3 ,^3)hl- 7 r + 2/(Ml - U 3 , V 2 )f(u 3 , V 3 ) In 



^l(M3 - M3J *" M 3 (^Mi - f 3 ) ' 

with = 1 — v i and = 1 — «j. The remaining functions can be deduced making use of the 
symmetry relations (J5.8j) . Here the notation was introduced for the function 

/(«,») = -—, (5-14) 
V v — u 

which we already encountered at one loop, see Eq. ()4.17j) . 

A few comments are in order. According to ()5.12|) . the third region in Fig. El produces a 
vanishing contribution to the three-particle kernel. This comes about as a result of cancellation 
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of the contribution of the diagram shown in Fig. |2] (a) against the one coming from a mirror 
symmetric diagram. 

It is known from the QCD calculations that, in general, two-loop evolution kernels involve 
special functions (dilogarithms Li 2 (x)). Still, the two-loop result for Y^iulv) is represented in 
terms of elementary functions only: the F— functions are given by a product of one-loop elements 
1)5.14)1 dressed by logarithmic dependence on the light-cone momenta. 

The function f(u, v) has a pole at u = v. In the expression for the one-loop evolution kernel, 
Eq. (j4.17j) . this pole is regularized by the '-(-'-prescription. To two loops, the F— functions involve 
the product of two /—functions and develop poles at Ui = V\ and u 3 = v 3 . It turns out that these 
poles cancel in the sum f)5.7j) in such a way that the evolution kernel Y^iulv) only has integrable 
singularities at u\ = v\ and u 3 = v 3 . It is easy to see that the two singularities correspond to 
the limit when the internal gluons in the diagram shown in Fig. El (a) have vanishing light-cone 
momenta, k + = 0. We remind that, in the light-cone axial gauge A + = 0, the gauge field 
propagator (j4.1|) is not well-defined for k + = and has to be supplemented by a particular 
prescription, Eqs. ()4.2|) and ()4.3j) . The fact that the two-loop kernel f)5.7j) is integrable at u\ = v\ 
and u 3 = v 3 implies that Y^2 3 (u\v) is not sensitive to the choice of the axial gauge prescription 
and, most importantly, it has finite moments fldu^Y^iulrfu^u^u 1 ^. 

We remind that the entire two-loop kernel ()5.2|) has to fulfill the polynomiality condition, that 
is, its moments ought to be polynomial in the v— variables, Eq. 1)3.22)1 . It turns out that the three- 



particle kernel v[ 2 3(w|t;) alone does not obey this condition. Its moments are given by rational 
functions of f 1,2.3 decorated by logarithms and dilogarithms depending on the ratio Vx/v 3 and 
V2/V3. As we will demonstrate in Sect. I5.2~3*j the "unwanted" terms violating the polynomiality 
of Y^siulv) cancel in the right-hand side of Eq. (J5.2j) against similar terms coming from the 
moments of two-particle irreducible kernels V^, . 

To separate the "unwanted" terms it proves convenient to introduce the double plus-distri- 
bution. It represents a natural generalization of (j4.15|) for three-particle irreducible kernels and 
is defined as 



1 



YY> 3 (u\v) 



++ 



1 

[du} 3 [tp(ui, u 3 ) - tp(v 1: u 3 ) - (p(m, v 3 ) + (p(v!,v 3 )] V^(u|«) , (5.15) 



where (p(ui,u 3 ) is a test function. Since the kernel Vj^t/wlv) lives on the simplices (|4.1fij) and 
(j4.18|) . as exhibited by the step-function structure, the constraint on the u— momentum fractions 
arising from the integration measure can be omitted, i.e., 6(1 — u\ — u 3 )N^2 3 {u\v) = Yi2 3 (u\v). 
One can verify then that the kernel [V^3(' u l' t ')]++ defined in this manner satisfies the polynomi- 
ality condition. By definition, 

r(i) 



Y$ 3 (u\v)} +5(ui-vi) f d Ul Y^ 3 (u\v) (5.16) 
J ++ Jo 

+ 5(u 3 - v 3 ) j du 3 Y$ 3 (u\v) - 5{ui - vi)5{u 3 - v 3 ) / [du] 3 Y^ 3 {u\v) , 



where the integrals entering the subtraction terms are well-defined and can be evaluated with 
the help of ()5.7|) (see Eqs. ()B.3j) - (jB.5j) in Appendix |Bj. We conclude that the polynomiality 
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Figure 4: Two-loop corrections to the two-particle evolution kernel involving only gluons in 
internal lines. The mirror symmetrical diagrams are not displayed. The gluon polarization 
insertion is displayed in Fig. El 




Figure 5: One- loop gluon polarization contributing as an insertion into the one- loop gluon 
exchange diagram. The gluon and fermion bubbles in Fig. EJ(f) are be denoted as (f) g i and (f)f er , 
respectively. 

is violated by the subtraction terms in the right-hand side of (j5.16|) . The first two subtraction 
terms preserve the momenta of a single particle and take the form of two-particle irreducible 
kernels. The last subtraction term has the structure of a single-particle contribution with the only 
difference that it involves the integral f[du] 3 ^U3( u \ v ) which is a rather complicated function 
of all three momentum fractions, f 1,2,3- Later in this section, we will demonstrate that the 
subtraction terms are compensated by similar polynomiality breaking terms coming from genuine 
two-particle evolution kernels. 

5.2 Two-particle contributions 

Let us now consider the two-particle irreducible contribution to the two-loop evolution kernel 
(15. 2|) . It is determined by Feynman diagrams shown in Figs. 0] and El In distinction with the 
three-particle kernel, the two-particle kernel is not universal and depends on the gauge theory 
under consideration. For the sake of simplicity we will distinguish between QCD-like diagrams, 
that is, those involving only gluons, and diagrams with scalars. The former are relevant for both 
QCD and SYM theories whereas the latter only appear in M = 2 and Af = 4 theories. Let us 
analyze the two sets of diagrams in turn. 

5.2.1 Gluon contribution 

The diagrams with gluon exchanges are shown in Fig. HJ As was already mentioned, similar 
diagrams define the two-loop contribution to the evolution kernel for the twist-two maximal 
helicity operators in QCD defined in (j2.18|) . The important difference between the three-particle 
operators, Eqs. (|2.2Uj) and (|2.19j) . and the two-particle operators, Eqs. (|2.18|) and (|2.17|) . is that in 
the latter case the pair of fermions is in the color-singlet state. However, this only affects the color 
factors accompanying the individual Feynman diagrams shown in Fig. HJ not their momentum 
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structure. Let us examine these color factors and relate to each other the two-particle contribution 
to the dilatation operators for the two- and three-particle operators. 

To start with, we present the result for the two-loop evolution kernel for the twist-two oper- 
ators (j2.18|) in QCD. In the momentum representation, it has the following form for V\ + v 2 — 
U\ + u 2 = 1 



^qcu(u 1 ,u 2 \v 1 ,v 2 ) 



9 2 C 



4vr 2 



rVS D («i.^i) + 



9 2 C F 
4tt 2 



V§cd(«i, «i) 



(5.17) 



where Cp = (N^ — 1)/(2N C ) is the Casimir operator in the fundamental representation of the 
SU(N C ). To one- loop order, the evolution kernel is determined by the universal expression ()4.17j) 



V&(u,v) 



u 8(v — u) 



u — ► u 

V — > V 



(5.11 



V v — u 

To two-loop order, the kernel Vqc D (w, t>) has been first calculated in Ref. [HZ] in the so-called 
forward limit, for v\ + v 2 = U\ + u 2 = 0. It was later reconstructed for u\ + u 2 ^ 0, making use 
of the conformal symmetry constraints, and yielded the result 



C F - hN, 



2 iV ^r(l) 



Pi 



C h 



C 



u — > u 



v — > V 



(5.19) 



where /3 QCDj o = llN c /3 — 2rif/3 is the one-loop beta-function in QCD, Eq. (|2.8J1 . and 



Wp\u,v) = — lnwlnw 



2v 

+ 6(v - u) 
Y^\u,v) = 9(v-u) 



9(v-u) 



(5.20) 



f(u,v) ( g-C(2) 



3 , u 1 . u v — u 
— In In — In 

4 v 2 v u 



1 u v-u' 
-/(w,f)ln-ln 

2 f f 



— Li 2 (w) — Li 2 (t> ) I + f(u, v)\nu\nv + 



u 

2v 



(5.21) 



f(u,v) ^Li 2 — ~~) + Li 2 (u) + ^ln 2 w - In it In' 

(li 2 (u)-Li 2 (l- 





u 







V^ 1} («,v) = 6{v-u) 



In 5 
-In — | 

4 w 12 



f(u,v) 



(5.22) 



Here £(n) is the Riemann zeta function, Li 2 (x) = — J* y ln(l — t) is the Euler dilogarithm and 
the /—function is defined in ([5. 14)1 . 

Comparing the QCD expression ()5.17|) with Eqs. (|3.2|) . (|4.10j) and (|4.17|) . we see that to 
one-loop order the evolution kernel Vi 2 for three-particle operators can be deduced from (J5.17j) 
through substitution of the color factors 



C f 



(5.23) 



for fermions in the fundamental and adjoint representations, respectively. Unfortunately, this 
simple property gets lost starting from two loops. To see this, one examines the color factors 
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corresponding to the Feynman diagrams shown in Fig.|Uin four different cases: twist-two operator 
(|2.18|) in QCD, its counter-part (|2.17|) in SYM theory, three-quark operator (|2.2()jl in QCD and 
three-gaugino (J2.f9|) operator in SYM. 
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We notice that for two-particle operators, going over from QCD to SYM amounts to replacing 
C F by the Casimir operator in the adjoint representation C F —>■ N c and substituting the number 
of quark flavors as nj- -» N c Af (see Eq. (I2.fij) ). Let us now turn to three-particle operators and 
substitute the color factors in the second row of the table according to (j5.2r>j) (we recall that the 
three-quark operators are well-defined for N c = 3 only). It is easy to see that this leads to the 
correct expressions for the color factors in the third and fifth rows except for the columns (d) 
and (e). Thus, in order to reconstruct the two-particle contribution to the evolution kernels for 
the three-quark and three-gaugino operators from the QCD expression (|5.19j) with the help of 
(I5.23|) . it suffices to supplement (|5.19j) with the contribution of two additional diagrams in QCD 
shown in Fig. |U(d) and (e). More precisely, one only needs the abelian part of their contribution, 
i.e., the one proportional to C\. In what follows, we denote the corresponding kernel as Vsv- A 
straightforward calculation leads to 

V$(«, v) = 0(v - u) f(u, v) ( - - In - - - In 2 - - - In - In {v - u)u (5.24) 

\ 4 v 4 v 2 v vv 

1 fu-v\ 1 fv-u\ 1 . 1 . 

+ -Li 2 — - -Li 2 + -Li 2 (u) - -Li 2 (v) - 





2 V v J 2 

, . V — U . U 1. 2 « 1, 

+ t(u, v) -In In In mumv 

1 2 v v 4 v 2 

1 T . / x 1 T . / x 1 T . fv — u\ 1 T . fu-v 
+-U 2 (u) - -U 2 (v) + ~U 2 [—)-^ 2 

u — > U ' 
V — > V 

We recall that ()5.17|) is defined for v\ + v 2 = ui + u 2 = 1 while for three-particle operator the two- 
particle evolution kernel, say , the momentum conservation implies that v 1 + v 2 = U\ + u 2 = 
1 — V3. To verify this condition, one can rescale the momentum u and v in (|5.17|) making use of 

(EE}. 

In this way, we obtain the following expression for the gauge field contribution to the two- 
particle kernel, > coming from the diagrams shown in Fig. 0] 

VgluoB^i^s \vi,v 2 ) = ■ (V^ + VW) : , : , 5.25 

8 Ml + U 2 \ J V Ul + U 2 Vi + f 2 / 
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Figure 6: Two-particle irreducible contribution involving scalars in internal lines. 



where Vgy is denned in (J5.24|) and Vqc D is given by the QCD expression (J5.19|) with the color 



factors substituted as 



QCD : 
SYM : 



C F 

c F 



N c /2, n f ^N c M 



(5.26) 



Notice that the — term, Eqs. 1)5.19)1 and (|5.21jl . only contributes in QCD and vanishes for 
gaugino in the adjoint representation. As we will show in Sect. El this leads to a dramatic differ- 
ence in integrability properties of the two-loop dilatation operator in QCD and SYM theories. 

The two-particle kernel Vgiuge does not depend on the renormalization scheme employed — it is 
the same within the regularization by means of DREG or DRED provided that the perturbation 
series is in terms of MS— coupling constant. Since it more natural to define the coupling constant 
in regularization via DRED in the DR— scheme, the two-particle dilatation operator acquires an 
extra term as explained below in Sect. 15 .2. HI 



5.2.2 Scalar contribution 

The scalars contribute to the two-particle evolution kernel in the M = 2 and M = 4 SYM theories. 
The corresponding Feynman diagrams are shown in Fig. El As before, the ladder diagram in Fig. 
El (e) vanishes due to the maximal R— charge of the three-fermion operator. The color factors 
of the diagrams in Fig. El (a), (b), (c) and (f) coincide with those shown in Fig. |2] (c), (d), (e) 
and (b), respectively, and can be read from the Table given above. In particular, the diagram in 
Fig. El (f) has a vanishing color factor and does not contribute. The remaining diagrams take the 
form of the one-loop graph (see Fig. ^) with the propagators and the vertices "dressed" by the 
scalar interaction. Therefore, it is not surprising that their contribution is proportional to the 
one-loop /—function ()5.14j) decorated by additional logarithms 

+0(v 2 - u 2 ) (in ^-fjf K v 2 ) I , (5.27) 

with n s = 2{M — 1) being the total number of scalars, Eq. ()2.6)) . The kernel ()5.27)) is the same 
in the MS— and DR— schemes. 

Notice that the scalar contribution ()5.27)) is similar to the Vp— term, Eq. ()5.22|) . In the 
gluon contribution, Eqs. (j5.19|) and (|5.25|) . this term is accompanied by the prefactor ~ /3 QC d,o = 
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lliV c /3— 2n//3, which coincides upon substitution ()5.26|) with the one-loop beta-function in SYM 
theory modulo the scalar contribution given by —n s /6. Combining together ()5.25j) and (|5.27jl . 
one finds that the logarithmically enhanced part of the scalar contribution (|5.27|) is accompanied 
by the factor +n s /Q. Writing this factor as —n s /6 + n s /3 and splitting the scalar contribution 
into two terms, one finds that the first term forces the coefficient in front to the Vp— term to be 
equal to the complete one-loop beta-function in the SYM. As we will show in Sect. 16.1.41 the 
second term is an artefact of the MS— renormalization scheme used in the evaluation of ()5.27j) . 

5.2.3 Two-particle kernel 

Let us combine together various two-loop contributions to the two-particle kernel , Eq- ()5.2|) . 
They come from three different sources: the gluon kernel f)5.25|) . scalar kernel ()5.27j) and two- 
particle kernel induced by the subtraction terms in (|5.16|) 

Vg = V« G >i, uah, v 2 ) + V^ la >i, UaK v 2 ) + V« K u 2 \ Vl , v 2 ) . (5.28) 

To identify the last term one substitutes the three-particle kernel ()5.16|) into (|5.2j) and selects the 
two-particle contribution proportional to 5{u 3 — v 3 ) = 8{u\ + u 2 — v x — v 2 ). In this way one gets 

V sub( M i' M 2|fi,f2) = S(ui +U2-V1- v 2 ) [W(u 1 ,u 2 \v l ,v 2 ) + W(u 2 ,ui\v 2 ,vi)} (5.29) 

where ^ 

W(ui,u 2 \vi,v 2 ) = du' 3 Y$ 3 (ui,l-u 1 -u' 3 ,u' 3 \v 1 ,v 2 ,v 3 ) (5.30) 





- u 2 1 




- v 2 J _ 



with u 2 = v\ + v 2 — u\ and v\ + v 2 + v 3 = 1 . The explicit form of V^ b (ui, u 2 \v i, v 2 ) can be found 
in Appendix iBl 

The two-particle kernel ()5.28|) . contrary to its three-particle counterpart, is not universal and 
depends on the gauge theory under consideration. Substituting ()5.25|) . ()5.27|) and ()5.29p into 
(I5.28J) one finds after some algebra 

Via = S(ux + u 2 - vx - v 2 ) 9(u 2 )9(u 2 - vx - v 2 )J r i(u\v) (5.31) 

5 

+9(ui)6(vi - ui)T 2 {u\v) + -9(u 2 )9(vi - u 2 )J 7 3(u\v) 

8 

where Uj = 1 — Uj and we used the results flB.8|) - (jB.9|) from Appendix El The {. . .}— term 
ensures symmetry of the two-particle kernel under interchange of two particles. The .T 7 — functions 
which enter as coefficients in front of different step-functions arise either from the subtraction 
terms alone (J-"i) or the genuine two-particle irreducible kernel (JF 3 ) or the sum of both {J 7 ^). The 
former being defined by the three-particle irreducible kernel is given by a universal expression 

.Fi(ti|v) = ^ r - \L{-u\, vx - v 2 , v 2 - u 2 ) - \L{-ux, vi - v 2 , vx) (5.32) 

{vx-v 2 ){vx-ux) 

+ \L(u 2 - v 2 , v 2 , vx - v 2 ) + \L(u 2 - v 2 , v 2 , vx) , 
with L(u, v, w) = ± In s [4 + In *j ] . The two remainine functions are defined as follows: 
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In SYM theory with M supercharges one has 



F<2{u\v) 



FoOulv) + 







i ( Pn,c 



N r . 



+ N 



ln^ + * 
V\ 3 



/K Vi)-(JV -l)/(«i, (5.33) 



where /3jv\o — N c (4 — A/") is the one-loop beta-functions in SYM, Eqs. (|2.9jl . and 

•FoMu) = \L(u 2 - v 2 , v 2 , vi - v 2 ) + \L(u 2 - v 2 , v 2 , vi) - \L(u 2 - v 2 , v l7 vi - v 2 ) 



(5.34) 
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with the /—function defined in (|5.14J1 . 
• In QCD the same functions are given by 
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(5.35) 
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«1 \ , «2 



2Li 2 1 - 



«2 

i; 2 



+ ln 2 



V-2 



Vi + v 2 



Vi + v 2 



In 



v 2 



In 

Vl +V 2 Vi+ v 2 
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where /3 QCDj0 = 11 — In/ is one-loop beta- functions in QCD, Eq. (|2.8|h and 



-^qcd = \ 

3 t>i(fi + w 2 J 



+ 2/(u 2 ,v 2 )ln 



u 2 



Vl 



In 

^1 + u 2 f 1 + w 2 
«i 



(5.36) 
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The following comments are in order. 

The step-function structure in front of T\ looks quite unusual and even appears not to con- 
tribute to the moments ()3.22|) which are defined with the measure [du] 3 possessing the support 
region {0 < Ui < l,$2i u i = 1}. If one uses the momentum-conserving condition from the 8- 
function accompanying T\ to express 6{u 2 — v% — v 2 ) as 9(— Ui) one would immediately come 
to this conclusion. This is not the case however since this contribution emerged from the sub- 
traction term in the definition of the double plus-distribution in Eq. (|5.16|) . As it is clear from 
there, while the Vi 2 3— kernel itself resides on the simplices (|4.16|) and (|4.18|) . the subtraction 
terms "live" in the entire domain of Ui. Their contribution to the two-particle subchannels was 
reconstructed from the subtracted expression by "integrating-in" corresponding 5— functions and 
assuming that the support of the resulting subtracted kernel coincides with the original unsub- 
tracted one. Therefore, the inconsistency which arises for the T\— term only is simply resolved 
by ignoring the spectral constraint on the total sum of momentum fractions after integrating out 
the momentum-conserving 5— function in the measure [du\ 3 . Namely, 



/ [rfM] 3 5(ui+u 2 -fi-f 2 )6 ) (M 2 )^(u 2 -fi-i; 2 )J'i(tt|v) = / du 2 Ti 

Jo J V1+V2 



(vi+v 2 -u 2 , u 2 , v 3 \vi, v 2 , v 3 ) , 

(5.37) 
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with v 3 = 1 — v i — v 2 on the right-hand side. Analogously, one treats other terms with the T\ 
structure. 

The expression (J5.31|) was obtained within the MS— renormalization scheme. Going over to 
the DR— scheme, one finds that contribution of all diagrams except the one containing gluon 
self-energy, Fig. 0] (f), remain the same. For the latter diagram, the difference between the 
two schemes arises at the level of 0(e) corrections to the gluon self-energy (see Fig. |3J). These 
corrections interfere with a double pole coming from the Feynman integral to produce a nontrivial 
contribution to the two-particle evolution kernel which is proportional to the one-loop evolution 
kernel. We will return to this issue in Sect. 15.4.11 

Comparing (|5.33|) and (|5.35|) one notices that the two-particle kernel in SYM theories has 
a simpler form and, most importantly, it does not involve the Li 2 — function. This property is 
nontrivial since Li 2 enters into the expression for gluon kernel (|5.25|) through the abelian contri- 
bution Vgy, Eq. (J5.24)) . A close examination reveals that in the expression for the two-particle 
kernel (|5.28|) the Li 2 — terms coming from gluon kernel cancel against similar terms generated by 
subtraction term ()5.29|) . Similar cancellation takes place inside the QCD expression ()5.35|) but 
in that case Li 2 — term also comes from the nonplanar kernel V^. In SYM theories this term 
does not contribute to the three-particle operators since it is accompanied by the color factor 
Cf — N c /2 which vanishes upon the substitution ()5.23|) . 

As was mentioned in Sect. |21 the conformal symmetry is broken in gauge theories with nonva- 
nishing beta-function. One observes that the kernels (J5.33|) and (J5.35j) contain terms proportional 
to the one-loop beta-function. In the light-like axial gauge, they come from the one-loop Feyn- 
man diagram in which the internal gluon gets "dressed" by a one-loop self-energy correction (see 
Fig. |H (f)). The corresponding terms are responsible for explicit breaking of conformal invari- 
ance of the evolution kernels at two loops. Other more subtle symmetry breaking effects will be 
addressed below in Sect. 16.1.41 

In distinction with the three-particle irreducible kernel (15 ,7j) . the two-particle kernels (J5.33|) 
and ()5.35|) are singular for Uj — > Vj. Carefully examining (|5.31|) in this limit, one finds that the 
two-particle kernel develops simple poles ~ — Vj) and, as a consequence, its moments 

are divergent. In the axial gauge with the principal value prescription, these poles are regularized 
according to (J4.3j) and lead to ~ In 5 contribution to the moments. These singular terms can 
be separated with the help of the plus-distribution as in Eq. f)4.15j) . For the pair-wise kernels 
— V^(«i, u 2 \vi, v 2 ) it reads 

[V (1) (wi, u 2 \v 1, f 2 )]+ = V (1) («i, u 2 \vi, v 2 )-5(ui-vi)5(u2-v 2 )j- duidu 2 Y ( - 1 \ui, u 2 \vi, v 2 ) , (5.38) 

where the integral in the right-hand side is evaluated using the principal value prescription. The 
two-particle kernel [V^ ]+ defined in this way has finite moments and satisfies the polynomiality 
condition. 

In the SYM theories with M supercharges, the M— dependence enters into the two-particle 
kernel ()5.31|) through the last two terms in the right-hand side of ()5.33|) . Taking into account that 
@N ,0 = (4 — M)N C one finds that the M— dependence only resides in the last term proportional 
to the one- loop /—function, Eq. (|5.14j) . Together with (J4.17j) this implies that the two-particle 
kernels in SYM theories with different number of supercharges are related to each other in a 
simple manner 
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Figure 7: Two-loop fermion wave functions renormalization in QCD and M = 1 SYM. 



where = V^(iti, u 2 \vi, t> 2 ) is the one-loop expression, Eq. (j4.17|) . This suggests that the 
two-loop dilatation operators for different M should be related to each other in a simple manner 
(see Eq. (jEEH) below). 

5.3 One-particle contributions 

Finally let us consider single-particle irreducible contributions V- to the dilatation operator 
(J5.2|) . They come from the fermion wave function renormalization and are defined by the Feynman 
diagrams shown in Figs. U\ and El 

The two-loop renormalization constants ^ in Eq. (|5.4|) depend on the renormalization 
scheme used for their evaluation. In addition, in the light-like axial gauge A + (x) = sup- 
plemented with the principal value prescription f)4.3|) . they depend on the ratio of the light-cone 
momentum carried by the fermion and the regularization parameter S 

( 5 - 4 °) 

The result of the calculation in the MS— scheme of ^ 1 \v/5) both in QCD [HH] and SYM theories 
is 

+ |_3f(2) + B{ (S)-^(l + |c(2) + ^]n£), (5.4!) 
d 1, (*+/«) = l+3C(3)+( 1 |-2C(2))lni- + !|_, ! ,^li +hl i_^ . (5.42) 

We recall that for three-quark operators in QCD the color factors take the values Cp = 4/3 and 
N c = 3. 

In Eq. (|5.2|) . the self-energy correction (|5.4U|) is accompanied by the delta-functions 5{ui — 
V\)8{u2 — v 2 ). Similar contribution also comes from the subtraction terms entering into the defi- 
nition the double-plus distribution for three-particle kernel (|5.16J) and from the plus distribution 
for two-particle kernel ()5.38|) . Combining them together one gets 

r« = i J [du] 3 Y^u\v)-^du 1 du 2 Y^(u 1 ,u 2 \v l ,v 2 )-^ [^\ Vl /5)+^\v 2 /5)]+... (5.43) 

where ellipsis denote the terms needed to restore the symmetry under cyclic permutations of 
particles (vi,v 2 ,v 3 ) — > (v 2 ,v 3 ,vi) — > (^3,^1,^2)- Two comments follow. First, notice that the 
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Figure 8: Two- loop Feynman graphs for wave function renormalization involving scalar fields 
propagating in internal lines. 



integrals of three- and two-particle irreducible kernels in the right-hand side of (|5.43jl contain 
a nontrivial dependence on momentum fractions fj. However, in the combination (|5.43J) this 
dependence entirely cancels as a manifestation of the restoration of the polynomiality condition. 
Second, the contributions of the two-particle kernel and self-energy corrections contain an explicit 
dependence on the axial gauge regularization parameter 5. The gauge invariance requires that 
the divergences ~ In 5 should cancel in their sum. Indeed, one can verify that stays finite 
for S — > 0. It reads in the MS scheme for the QCD case with N c = 3 



QCD 



27 13 

16 32 QCD '° ' 



(5.44) 



and for the SYM case 



1 



1 



-n: 



(5.45) 
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with n s = 2(7V— 1). By the construction, 3T^ defines the two-loop correction to the anomalous 
dimension of the local three-fermion operators, Eqs. ()3.13|) and ()3.15|) . in the MS— scheme. The 
one-loop correction is given by ()4.12|) . 



5.4 Two-loop result 

Summarizing the calculation, the two-loop correction to the dilatation operator ()5.2|) reads 
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-3T^6(u 1 -v 1 )S(u 2 -v 2 ), 

(5.46) 

where each term in the right-hand side has finite moments and satisfies the polynomiality con- 
dition. Here the three-particle kernel V^ 3 is defined in ()5.7|) . the two-particle kernel vj^ is 
given by (J5.31|) whereas the remaining kernels can be obtained from these two with the help of 
(|5.3j) . The '++'- and '-l-'-distributions were introduced in (|5.16|) and (|5.38|) . respectively. The 
single-particle contribution is specified in ()5.44|) and ()5.45|) . 



5.4.1 Scheme dependence 

The above mentioned results were obtained within the MS— renormalization scheme. As was 
already mentioned, the two-particle and single-particle contributions have different form in the 
MS— and DR— schemes. In QCD both schemes are equally suitable while in the SYM theory the 
DR— scheme is advantageous as it preserves supersymmetry. Therefore, as a final step we have to 
transform ()5.46|) to the supersymmetry preserving dimensional reduction scheme. This scheme 
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transformation is achieved via a finite renormalization of the nonlocal light-cone operators ()2.15jl 
renormalized according to the conventional MS scheme. For nonlocal operator of arbitrary length 
L it reads 

Obs(zi, ...,z l ) = w(g 2 ) Om(z u ...,z L ), (5.47) 

where w(g 2 ) is, in general, a complicated integral operator acting on z— coordinates. Substituting 
this relation into the evolution equation ()2.13|) one finds that the evolution kernels in the two 
schemes are related as 

M m (g 2 ) = w(g 2 ) M m (g 2 ) «T V) - /% 2 )^P «T V) . (5.48) 

Notice that in the conformal limit, for (3{g 2 ) = 0, the two evolution kernels are related to each 
other through similarity transformation and, therefore, have the same eigenvalues. The latter 
determine the anomalous dimension of Wilson operators in the underlying gauge theory. 

Equation ()5.46|) defines the two-loop correction to H^^g 2 ) in the momentum representation. 
In order to transform it into the DR— scheme with the help of (I5.48|) . it suffices to determine 
w(g 2 ) to one loop. To this accuracy, w arises due to the fact that the one-loop diagrams shown in 
Fig. Q have different finite parts when renormalized in the DR— and MS— schemes. It is easy to 
see that for maximal-helicity operators, w only receives contribution from the gluon self-energy, 
Fig. E (b), while the diagrams (a) and (c) are the same in both schemes to the ~ e° accuracy. 
Therefore, w is given by the unit operator times a factor. The latter can be easily computed 
with the result 

Yet in the above scheme transformation, the gauge coupling is still defined in the MS scheme. 
Thus, one has to further renormalize the coupling according to the DRED conventions, this yields 
its finite renormalization to two-loop accuracy jlU] 

9k = 9k(l-^&] • (5-50) 



6 8tt\ 

Thus, for the evolution kernel in the MS— scheme given by 

M m = AmsH (0) + A^H« + 0(\ 3 ) , (5.51) 
the same kernel in the DR— scheme takes the form 



or equivalently 



im = (l - ^p) M (0) + A^H«] - ^Als + 0(A 3 ) , (5.52) 

Hdr = A^eW + A^ (W) - V) - ^ + G(A 3 ) . (5.53) 

According to (J5.52|) . the two-loop dilatation operators in the two schemes differ by an overall 
normalization only and, therefore, they share the same symmetry properties (if any). This 
property is a unique feature of the dilatation operator in the maximal helicity sector and it holds 
in SYM theory independently of the number of supercharges M . 
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At L = 3, one finds from ()5.53j) that the two-loop evolution kernel for three-particle operators 
in SYM theory is given in the DR— scheme in the momentum representation by 

y« = yd) _ ly(0) _ 3^ _ ^ (5 _ 54) 

6 4 

with V(°) and defined in Eqs. ()4.13|) and ()5.46|) . respectively, while 11 is the identity operator 
11 = 5(ui - v\)8(u2 - v 2 ). 



5.4.2 Limiting cases 

The two-loop evolution kernel V*- 1 -* is given by a rather complicated expression. There are however 
two limiting cases in which its form simplifies significantly. Let us examine for Uj — ► Vj 
(with j = 1,2,3). As was explained in Sect. 14.21 this corresponds to the limit when gluons 
exchanged between the fermions in the Feynman diagrams shown in Fig. 0] carry vanishing light- 
cone momenta Uj — Vj — > 0. We argued in Sect. 15. 231 that for Uj — > Vj the dominant contribution 
to comes from the two-particle evolution kernels [V^.]+ which develop poles ~ l/(uj — Vj). It 
is straightforward to verify using ()5.31|) that these poles are generated by terms in the expression 
for Y^jj involving the one-loop /—functions, Eq. (|5.14jh 

y(i) S ^ + u<2 _ Vi _ v ^ _ Ul )f(u^ Vl ) + e(v 2 - u 2 )f(u 2 , v 2 )} . (5.55) 

One recognizes the expression in the right-hand side as the one- loop kernel defined in f)4.17)l . 
Substituting ()5.55|) into ()5.46|) one finds that for Uj — v j — > the full two-loop evolution kernel 
is proportional to the one-loop contribution, Eq. (|4.13jh 

Y(u\v) = XY {0 \u\v) + X 2 Y w (u\v) + 0(X 3 ) U] = 3 ^T cusp (X) ■ Y (0) (u\v) + . . . , (5.56) 

where ellipsis denote terms subleading as Uj—Vj — > 0. One can show that ()5.56j) holds to all orders 
of perturbation theory |H] with r cusp (A) being the universal cusp anomalous dimension [32] 
related to a cusp anomaly of Wilson loops |4l5] . The two-loop expression for r cusp (A) can be 
obtained from (|5.31|) . In SYM theory, in the DR— scheme one finds [4*4*] 

-21 
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2(4- AT) 



7T 



(5.57) 



with the coupling constant A = g^-N c / (8tt 2 ) . 

The second limiting case corresponds to the selection of terms inside Y(u\v) proportional to 
the one-loop beta functions. These terms are gauge invariant and originate in the axial gauge 
from a special subset of Feynman diagrams involving gluon propagator dressed by one-loop 
corrections (see Fig. HJ (f)). In QCD one can automatically select the required terms by going 
over to the limit of large /3o, that is g 2 /3o = fixed for /3q — > oo. Since /?o = 11 — 2ri//3, the large 
[3q— limit can be formally achieved in QCD by continuing Y(u\v) to the number of quark flavors 
71/ — > — oo. Similar procedure can not be performed in the SYM theory since the number of 
gaugino equals the number of supercharges, Eq. (|2.6|) . According to (|5.31|) . the /3 — term only 
comes from two-particle kernels and has the same form in QCD and SYM theories 



V 12 = TTT^K +U2-V!- v 2 ) 



^(»jWn-»i)(ln^ + | > )/(«i,« l ) + ( , * 1 S Ua 

vi 3/ \v 1 ^v 2 



+ ... 

(5.5J 
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where Cr = N c in SYM theories and Cr = 4/3 in QCD. A unique feature of fio enhanced terms 
is that they can be resummed to all loops. To n th order of perturbation theory, these terms 
take the form X(X(3 ) n and come from the one-loop Feynman diagrams in which the bare gluon 
propagator is dressed by "renormalon" chain involving n one-loop polarization operators (see 
Fig. Ej). These diagrams contribute to the two-particle kernels and lead to the following all-loop 
expression for the evolution kernel in the large (3$— limit 



V(u\v] 



V12 + V 23 + V 
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where A = g 2 Cn/ (87r 2 ) and the notation was introduced for the function 



<p(x) 
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(5.61) 

in powers of A are in 

agreement with ()5.58j) . Comparing ()5.6()j) with the one-loop expression for the evolution kernel 
()4.17|) one notices that higher order corrections "renormalize" the conformal spin j of the fermion 
fields and induce the additional normalization factor 



It is easy to verify that the first two terms of the expansion of 



This relation holds both in QCD and SYM theories. 
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5.4.3 AT— dependence 

To one-loop order, the kernel V(°) depends on J\f through the normalization constant T^ym; 
Eq. (J4.12j) . To two-loop order, in Eq. (J5.46|) . the M— dependence resides in the two-particle 



kernels and the normalization constant £ 
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Taking into account ()5.39j) . one finds 
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where the subscript refers to the number of supercharges. This leads to the following remarkable 
relation between the two-loop dilatation operators in SYM theories with different number of 
supercharges 

H Ar (A) =[1-{M - 1)A] • EW = i(A) + const + C(A 3 ) . (5.64) 

We conclude that to two-loop accuracy the spectrum of anomalous dimensions in the sector of 
maximal helicity operators looks alike in all SYM theories including the maximally supersym- 
metric M = 4 theory. If the dilatation operators in these theories have different properties, the 
difference could only appear starting from three loops. Eq. ()5.64|) is one of the main results of 
this paper. 
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6 Eigenspectrum of the dilatation operator 



In the previous section, we have calculated the two-loop dilatation operator in the momentum 
representation. Let us now turn to the analysis of its symmetry properties. Our strategy will 
be to solve the spectral problem for the two-loop evolution kernel Y(u\v), Eq. 1)3.25)1 . and, then, 
look for the hidden symmetry of its spectrum. 

As we already observed in Sect. 13.1) the dilatation operator is invariant under translations and 
dilatations along the light-ray, Eqs. (j3.12|) and 1)3.11)1 . respectively. The latter are the light-cone 
projections of the Poincare transformations in the underlying gauge theory. In the momentum 
representation, this symmetry leads to scaling property of the evolution kernel 1)3.21)1 which in 
its turn allows one to identify its eigenstates as homogeneous polynomials, Eq. 1)3.28)1 . It is 
well-known that classical Yang-Mills theory is invariant under an even larger set of space-time 
transformations — the £0(4, 2) conformal group. The conformal symmetry imposes additional 
constraints both on the properties of the evolution kernels and their eigenspectra. In particular, 
it allows one to resolve the mixing between the Wilson operators involving total derivatives 

On quantum level, the conformal symmetry becomes anomalous for /3(g 2 ) ^ 0. As we will 
explain later in this section, the conformal anomaly affects the dilatation operator starting from 
two loops only and, therefore, the one-loop dilatation operator, Eqs. ()4.7jl and (j4.13)l . inherits the 
conformal symmetry of the classical theory. We will demonstrate that the conformal symmetry 
allows one to diagonalize the two-particle evolution kernels in the coordinate and momentum 
representations, Eqs. 1)4.11)1 and ()4.17j) . respectively, but it is not sufficient to diagonalize the 
three-particle dilatation operators. 

Before we turn to the three-particle operators, let us examine the spectrum of the dilatation 
operators for the twist-two operators, Eqs. ()2.18)) and ()2.17jl . The reason for doing this is twofold. 
Firstly, the evolution kernel for twist-two operators contributes to the two-particle kernel V12 in 
the expression for the evolution kernel (J5.46jl . Secondly, this offers a convenient framework for 
discussing the conformal symmetry breaking. 

6.1 Twist-two operators 

The generating function for twist-two operators in QCD and SYM theories is defined in ()2.18jl 
and ()2.17jl . respectively. Similar to ()3.27j) . the twist-two operators are uniquely defined by poly- 
nomials P(ui, u 2 ) which are eigenstates of the evolution kernel in the momentum representation, 
V(u%, u 2 \vi, v 2 ). In QCD the two-loop evolution kernel Y QCD (ui, u 2 \vi, v 2 ) was defined in (15.17)1 . 
In SYM theories the evolution kernel Y SYM (ui, u 2 \v 1 , v 2 ) can be obtained from the QCD expres- 
sion by properly adjusting the color factors and adding the contribution of scalars determined in 
Sect. EH 

6.1.1 Eigenspectrum at one-loop 

To one-loop order, the evolution kernel for twist-two maximal-helicity operators in QCD and 
SYM are given by the same universal expression ()5.18jl 

VS- 2 («I«) = 2 [V (0) («i, «2K«2)] + - 2T^S( Ul - v x )8{u 2 - v 2 ) , (6.1) 
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with I^ ) defined in (I4.12J) and being the two-particle evolution kernel, Eq. ()4.14j) . Solving 
the spectral problem for this operator 

[du} 2 [N^\u l) u 2 \v l) v 2 )] + P^\u l) u 2 ) = - T (n)Pf { Vl ,v 2 ) , (6.2) 

with [du] 2 = duidu 2 5{u\ + u 2 — v i — v 2 ), it is straightforward to verify that its eigenstates are 
given by Gegenbauer polynomials j2Hl EZ| 

Ptf\ Ul ,u 2 ) = (m+u^cT ( V ^\ , (6.3) 

\U\ + u 2 J 

with n and I nonnegative integers. The corresponding eigenvalues do not depend on i and are 
given by the Euler digamma function 

7(71) = 2ty>(n + 2) - V(2)] . (6.4) 

The eigenstates of the two-particle evolution kernel, Eq. ()6.3|) . define the local twist-two operators 
O nl (0) with I > n 

O nl (0) = Pi l - n \d u d 2 )O(z u z 2 )\ zi=z2=0 , (6.5) 

where 0(21,22) is given by (j2.18|) and (j2.17|) . As follows from ()6.3|) . the operators O n i(0) with 
/ > n can be obtained from the one with I = n by applying the total derivatives, O n i(0) = 
d l + n O nn (0). 

The operators ()6.5|) have an autonomous scale dependence to one-loop order 



d 



,2 



»fO nl (0) = -|^ 7 (°)^(0) + 0(g 4 ) , (6.6) 
with the anomalous dimension given in QCD and in the SYM theory with M supercharges by 

7i 0) L. = C F [A^{n + 2) - 40(2) + 1] , (6.7) 



I QCD 

^ 0) | SY m = N c [4^(" + 2) - 4^(2) + A^] . (6 A 



Although these expressions define the anomalous dimensions for integer n > 0, they can be 
analytically continued to the complex n— plane. In particular, for n = — 1 one finds in the SYM 
case that the anomalous dimension coincides with the one-loop beta-function, Eq. (|2.9j) . 

7S-1 = -fafl = N c [N - 4] (6.9) 

As we will argue in Sect. 16.1.51 this leading order property is not accidental and can be extended 
to all loops. 



6.1.2 Conformal symmetry 

As was already mentioned, the explicit form of the one-loop eigenstates ()6.3j) is uniquely deter- 
mined by the conformal symmetry. We remind that the nonlocal light-cone operators 0(21,22) 
are built from elementary fields X{zn iJ ) located on the light ray, defined by the null vector n^. 
For such operators the full SO(4, 2) conformal group consisting of Lorentz rotations M M „, trans- 
lations P M , dilatation D and special conformal boosts K M is effectively reduced to its "collinear" 
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SL(2) subgroup which acts nontrivially on functions of the light-cone coordinate z — x ■ n*, 
defined by means of a tangent vector n* to the light cone such that n • n* — 1. The generators 
of the collinear subgroup are 

L - = -iP+, L + = |K_, L° = ^(M_ + + ©), (6.10) 

with the standard commutation relations 

[L°, it] = ±h ± , [L+ L - ] = L° . (6.11) 

These quantum field operators can be expressed in terms of the energy-momentum tensor in the 
underlying gauge theory. Still, for their action on the elementary field operators X(z k ) = X{z k n ll ) 
they can be represented by differential operators acting on the light-cone coordinate 

\L~,X(z k )] = L-X(z k ) = -d Zk X(z k ), 

[h+,X(z k )\ = LlX(z k ) = (z 2 k d Zk + 2z k j)X(z k ), (6.12) 
[h°,X(z k )} ee L° k X(z k ) = (z k d Zk +j)X(z k ), 

where the notation was introduced for the differential operators L k (with a = ±, 0) acting on 
the light-cone coordinate z k . Here j is the conformal spin of the field constructed from its scaling 
dimension d and the projection of the spin onto the light-cone s 

j = l(d + s). (6.13) 

For gaugino/quark operators defined in (|2.14jl and ()2.16|) . these parameters admit the values at 
the classical level d = 3/2 and s = 1/2 leading to j = 1. In quantum theory, the scaling dimension 
d acquires an anomalous contribution while s is protected to all orders of the perturbation theory. 

The field X(z k ) defines a representation of the SX(2;R) group which we shall denote as 
Vj. The representation space Vj possesses the highest weight X(0), such that L + X(0) = 
and L°X(0) = jX(0) while the descendant states are given by its derivative (L~) fc X(0) = 
(— l) fc <9^X(0). The nonlocal operators Q(zi,z 2 ) belong to the tensor product of two SX(2;R) 
modules which can be decomposed over the irreducible components as 

H®V, = ^2v 2j+n . (6.14) 

n>0 

The representation space of V^+n is spanned by local twist-two Wilson operators O n i(0). Among 
these operators one can distinguish the highest weight / = n satisfying the relations ^Hl HHj 

[L + , O nn (0)} = , [L°, O nn (0)] = (2j + n)O n0 (0) (6.15) 

and its descendants with / > n 

O nl (0) = [[O n „(0),L-],L-]...L-]] = (d + ) l ~ n O nn (0) . (6.16) 

The operators O n i(0) carry the conformal spin 2j + n and have the canonical dimension I + 3. 

It is straightforward to verify that the twist-two operators defined in ()6.5|) and ()6.3|) satisfy 
the relations ()6.15|) and (|6.16J) and, therefore, they belong to the SL(2; R) module V2j+ n - Then, 
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the conformal invariance implies that the Wilson operators belonging to different SL(2) modules 
do not mix with each other and, in addition, the operators belonging to the same module have 
the same anomalous dimension. 

Let us define the representation of the conformal group generators for the nonlocal light-cone 
operators 0(21,22) similar to ()6.12|) . Evaluating \L a , Q(zi, z 2 )) (with a — ±, 0) and making use 
of (j6.12|) . it is easy to see that to the lowest order of perturbation theory the generators h a are 
given by the sum of differential operators Lf + L 2 defined in (16.12)1 . In two of these operators, L~ 
and L°, we immediately recognize (up to an additive constant) the differential operators which 
were found to commute with the dilatation operator, see Eqs. ()3.11|) and (J3.12J) . The conformal 
symmetry implies that the one-loop dilatation operator also commutes with the conformal boost 
L + leading to 

[L? + ^,Hg] = 0, (6.17) 

with a — ±, 0. Notice that in distinction with (j3.11j) and (j3.12|) this property is lost starting from 
two loops due to the conformal anomaly. We remind that the evolution kernel Y^(u%, U2\vi, v 2 ) 
defines the operator M.f 2 , Eq. ()4.11|) . in the momentum representation. 

It follows from f)6. 17j) that the two-particle evolution kernel M.f 2 is a function of the quadratic 
Casimir operators defined on the tensor product ()6.14|) 

L\ 2 = (L° 2 ) 2 + i (L+ L7 2 + L^L+ ) (6.18) 
= -( Zl - z 2 ) 2 ^d Zl d z M - z 2 ) 2 ^ = J 12 (J 12 - 1) 

with L" 2 = (Li + L2) a and J\ 2 being the two-particle conformal spin. In the coordinate repre- 
sentation, the SL{2) module V 2 j +n is spanned by homogeneous polynomials 

^\z u z 2 ) = ( Zl -z 2 )\ ^\z u z 2 ) = (Ll 2 y^\z u z 2 ). (6.19) 

with the lowest weight, L^^r? = 0, and \&n its descendant. These polynomials are in the 
one-to-one correspondence with similar polynomials in the momentum representation, Eq. ()6.3|) . 
and satisfy the orthogonality condition similar to ()3.26|) . The eigenstates f)6. 19j) diagonalize the 
Casimir operator (J6.18|) and, as a consequence, the Hamiltonian (J4.ll)) has the same eigenvalues 
on these eigenfunctions independent of £, 

Hg.*W(zi,Z2) = 2[i;(n + 2j)-ij(2j)}^(z 1 ,z 2 ), (6.20) 
L\ 2 ^\z x ,z 2 ) = {n + 2j)(n + 2j-l)¥S{z 1 ,z 2 ) 

with ip(x) = dlnT(x)/dx. This relation is a counter-part of ()6.2|) in the momentum representa- 
tion. Together with J 12 • = (n + 2j)^n it allows one to write down the two-particle kernel 
in the SL{2) invariant form 

M^ = 2[^(J 12 )-^(2j)}. (6.21) 

Here j — 1 is the conformal spin of the quark/gaugino fields to the lowest order of perturbation 
theory. 

6.1.3 Conformal symmetry breaking 

Let us now examine the two-loop corrections to the dilatation operator for the twist-two oper- 
ators. As was explained in Sect. 15.4.11 starting from two loops the dilatation operator depends 
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on the renormalization scheme. Throughout the paper we employ the MS— scheme and its spin- 
off the DR— scheme. Another subtle issue emerging at two loops is the breaking of conformal 
symmetry in the dilatation operator within these two schemes even for (3(g 2 ) = 0. 

It is a matter of traditional folklore that in a generic gauge theory the conformal symmetry 
is broken if the beta-function is nonvanishing. This is a direct consequence of the dimensional 
transmutation phenomenon which generates an intrinsic mass scale in gauge theory modifying 
the scaling behavior of correlation functions. The nonzero beta-function induces an anomaly 
in the trace of the energy-momentum tensor. The latter is then generates in turn anomalous 
dimensions of composite operators when their product is renormalized in the conformal Ward 
identities, as we will demonstrate below. However, even in cases when the four-dimensional 
beta function is zero to all orders of perturbation theory, there is yet another source of the 
conformal symmetry breaking related to the choice of the regularization procedure. We remind 
that in the dimensional regularization with D = 4 — 2e the beta-function acquires an additional 
contribution (3 £ (g 2 ) ~ — 2e, Eq. (14. 6|) . and, therefore, the conformal symmetry is violated for 
e 7^ 0. According to ()4.4|) and (|4.5J1 . the dilatation operator is related to the residue of a simple 
pole 1/e in the expression for renormalized light-cone operator O(zj). Subtracting divergences 
and sending e — > afterwards, one generates symmetry breaking contribution to the dilatation 
operator coming from terms ~ /3 E (g 2 )/e. This source of the symmetry breaking is a peculiar 
feature of dimensional regularization rather than intrinsic property of the dilatation operator. In 
other words, in gauge theories with vanishing beta-function the conformal symmetry breaking 
terms can be removed by performing a scheme transformation of the dilatation operator and 
by going over to the so-called conformal scheme. Obviously, this transformation does not affect 
the eigenvalues of the dilatation operator but it does change the form of the corresponding 
eigenstates. 

The conformal operators O n /(0) define the basis of twist-two operators in gauge theory. To 
one-loop order these operators have an autonomous scale dependence but they mix starting from 
two loops 



and its form to higher orders is constrained by the Lorentz symmetry. To see this, one considers 
the operator O nn {zn lJ ) having the canonical dimension n + 3 and carrying the maximal possible 
conformal spin. It could mix with operators O mn {zn^) of the same canonical dimension but 
smaller conformal spin m < n. Notice that these operators necessarily involve (n — m) total 
derivatives. Then, expanding the operators around the origin z = 0, one finds that O n i(0) = 
d l + n O nn (0) could only mix with the operators Oki(0) with k < n. This means that the matrix 
of anomalous dimensions jnkig 2 ) has a triangular form, = for n < k, or explicitly 




(6.22) 




(6.23) 








\ 



7 



(i) 



(6.24) 
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with the diagonal entries jn = Inn and non-diagonal elements accompanied by the step function 
Onk — {^-i n > k;0,n < k}. The eigenvalues of the mixing matrix ()6.24|) are only determined by 
the diagonal matrix elements jn while the off-diagonal matrix elements only affect the form of 
the corresponding eigenstates. 

A convenient approach to study manifestation of the conformal symmetry breaking in the 
mixing matrix (|fi.24j) is through the conformal Ward identities. Below we briefly review this 
formalism, for a detailed discussion see j3U] and reviews |4*m i32j. Let us examine the variation of 
the correlation function (O n i(0) X(zi) . . .X{z£)) under transformations generated by the SL{2) 
generator L° 



U + Zkd k ) (O nl (0) X( Zl ) . ..X(z L )) = -(O nl (0) 5 (X( Zl ) . ..X{z L ))) (6.25) 

k 

= (5 O nl {0)X{ Zl ) . ..X(z L )) + (i8 S YM O nl {0)X{ Zl ) . ..X{z L )) , 

with 5 X(z k ) = — [L°, X(zk)] defined in (|6.12|) . The last term in the right-hand side is due to non- 
vanishing variation of the regularized Yang-Mills action. It is generated by a nonvanishing trace 
of the energy- momentum tensor in the D-dimensional gauge theory 8oS YM = —i f d D z O ^(z) . 
The product of the conformal operator and the trace anomaly requires an additional renormal- 
ization. This produces an anomalous contribution proportional to the anomalous dimension of 
the conformal operator O n i 



O nl (0)8 S, 



-7nk(g 2 )O kl (0) + . . . 



(6.26) 



where ellipsis denote "regular" terms. Here and in what follows the summation over the repeated 
index < k < n is implied. Combining together 1)6.25)1 and ()6.26jl one obtains that the conformal 
operator is transformed as 



[0,(0), L°] = ~ [(i + Aj)5 nk + lnk {g 2 )] O w (0) - % - 



d*zO uu {z) 



+ ... 



(6.27) 



where j is the conformal spin of the constituent field and the subscript 'R' indicates that the 
product of two operators is renormalized. In a similar manner, one can examine the variation of 
correlation functions under transformations generated by the SL(2) generator L + . They lead to 
the following commutation relations 



[O nl (0),L + ] = l - [a nk (l)+j c nk (l;g 2 )] O kl ^ - % - 



O nl (0) j (fz2z-0^{z) 



+ 



(6.28) 



where a nk {l) = a(n, l)5 nk = 2{n — Z)(n+/+4j — l)5 nk is a diagonal matrix and the matrix 7^ fc (/; g 2 ) 
is the so-called special conformal anomaly. This matrix has a well-defined perturbative expansion 
in powers of the gauge coupling and similar to the scale anomaly, Eq. ()6.26p . it arises from the 
renormalization of the product of the conformal operator and the variation of the Yang-Mills 
action, 

O nl (0)5+S YM = 1-^(1; 9 2 )Oh_i + . . . . (6.29) 



The one- loop expression for ^ c nk has been found in Ref. |STi] . 
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We remind that the SL(2) generators satisfy the commutation relation [L°,L + ] = L + . Eval- 
uating the commutators with O n i(0) on both sides of this relation and making use of (J6.27)) and 
(J6.28|) one obtains the consistency condition 

[7(9V + f(3 2 )+*>]=0, (6.30) 

which relates to each other the two anomaly matrices. Here, the a— matrix was defined in 
1)6.28)1 while the matrix b possesses the following off-diagonal elements b nk = —2(2k + 3)8 nk . In 
distinction with -y n k(g 2 ), the matrix elements 7^(/;5' 2 ) depend on the scaling dimension of the 
conformal operator /. This dependence can be fixed from the commutation relations [L°,ILr] = 
Lr leading to 

ri k (l + l; g 2 ) - i c nk (l; g 2 ) = -2 lnk {g 2 ) ■ (6.31) 

Since the one-loop correction to the matrix of anomalous dimension is diagonal, Eq. ()6.23jl . 
the relation ()6.3())1 is satisfied automatically to the lowest order of perturbation theory. Beyond 
leading order one finds that the special conformal anomaly has off-diagonal elements, 7^(/; g 2 ) ^ 
for n > k, and as a consequence the anomalous dimension matrix ^ nk {g 2 ) also acquires nonzero 
off-diagonal elements, even in theory with vanishing four-dimensional beta function. The reason 
for this apparently counterintuitive result is due to the fact that the beta-function is different from 
zero in D = 4 — 2e dimensions, the anomaly propagates into the finite result for the off-diagonal 
matrix elements of ■y nk (g 2 )- 

Expanding the left-hand side of ()6.30)) in powers of g 2 one observes that ()6.30)) mixes orders 
of perturbative expansion. Obviously, this relation does not provide any information on the 
diagonal elements of the anomalous dimension matrix 7„„ but it does constrain the off-diagonal 
elements jnk with n > k. In two-loop approximation, taking into account the explicit form of 
the one-loop special conformal anomaly matrix 7^(/; g 2 ), one finds the off-diagonal elements of 
the anomalous dimension matrix as [HI] 

n>k = ~ (7i 0) - l { k 0) ) {dnk{lt ] + A,) + <**}, (6-32) 

where the coefficients d nk are defined in terms of the derivative of the Gegenbauer polynomials 
over its index as 



^ = ~ 2 S dnkCriX) ' ^ = - (»-*)(* + * + 3) ' <6 ' 33) 



while the explicit form of g n k— coefficients is not relevant for our purposes and can be found in 
|5Uj . As we show in the next section, Eq. ()6.32|) allows one to restore the contribution of scalars 
to the two-loop evolution kernels in SYM without actual loop calculations. 

As we just explained, the conformal symmetry is broken in dimensional regularization even 
when the four-dimensional /3-function vanishes identically. In the latter case, the conformal sym- 
metry can be restored however by performing a finite renormalization of the conformal operators, 
Eq. O 

O nl (0) = B nk (g 2 )d kl (0) . (6.34) 

The £>-matrix absorbs symmetry breaking effects and rotates the conformal operators to a new 
basis, in which the operators O k i(0) do not mix and have an autonomous scale dependence (for 

P(g 2 ) = o) 

^d nl (o) = -in(g 2 )d nl (o) ■ (6.35) 
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The anomalous dimensions 7 n (g 2 ) coincide with the diagonal matrix elements of the mixing 
matrix (j6.23jl . 7 n = / y n n{g 2 )- The transformation (j6.34j) defines a new scheme known as the 
conformal subtraction scheme, or simply CS. By construction, the rotation matrix B diagonalizes 
the conformal anomaly matrix 

(B- 1 j(g 2 )B) nk = ln (g 2 )5 nk . (6.36) 

Wherefrom one can find that B~l = 6 nk + 7^(^ g 2 )/a(n, k). Then, it follows from ()6.30|) that, 
in the conformal limit f3(g 2 ) = 0, the B— matrix also diagonalizes the special conformal anomaly 
matrix 

(£T> + 7 V)] ®) nk = 2(n-l)(n + l + Aj-l + ln {g 2 ))8 nk . (6.37) 

In this limit, the energy-momentum tensor has vanishing trace, 9^ = 0, and, as a consequence, 
Eqs. (16.27)) and (|6.28|) simplify significantly. Combining them with (|6.34|) . (|6.36|) and (|6.37|) we 
conclude that the operators O n i(0) (with I > n) define a representation of the collinear SL(2) 
group. The highest weight of this representation, O nn (0), satisfies the relations 

[L+, 6 nn (0)} = , [L°, O nn (0)} = (n + 2j + i n (g 2 ))6 nn (0) , (6-38) 

which should be compared with (|6.15|) . One observes that, in the conformal limit, higher order 
corrections renormalize the conformal spin of the SL(2) module, to which the conformal operators 
O n i(0) belong to, by an amount proportional to their anomalous dimension. This property alone 
fixes the explicit form of the operators O n i(0). As before, the conformal symmetry constrains 
the form of twist-two operators but it does not allow us to determine the explicit expressions for 
their anomalous dimension. 

6.1.4 Contribution of scalars from one- loop special conformal anomaly 

Let us give a brief demonstration of the application of the conformal constraint (|6.32j) at two- 
loop order and how it compares with the conventional diagonalization procedure of the evolution 
kernels deduced in the previous sections. In QCD, the twist-two two-loop evolution kernel in the 
momentum representation was given in Eq. (j5.17J) . In SYM case, the two- loop evolution kernel 
receives an additional contribution from scalars determined by the Feynman diagrams shown in 
Fig. IH1 As was already explained in Sect. 15.2) the scalar contribution to the twist-two evolution 
kernel coincides, modulo the overall color factor, with the two-particle kernel 1)5.27), contributing 
to the three-fermion evolution kernel. 

In the basis of conformal operators O n i(0), Eqs. ()6.5j) and (j6.3j) . the scalar kernel (j5.27j) is 
represented by the mixing matrix 7^ s °. Similar to ()6.2|) . it can be found as 

/ [*.],2V£ 1 .(« 1 ,« a h,«)<r (^) = "E (iST/JV?) Of (^) , (6.39) 

with the measure [du\2 defined after Eq. ()6.2j) . Here, the factor 2 in the left-hand side accounts 
for the difference in color factors of singlet and octet kernels, cf. Table in Sect. 15.2.11 As expected, 
the mixing matrix has a triangular form, 7^ hC = for k > n. Making use of (|6.2j) . it is easy to see 
that the off-diagonal matrix elements 7^ sc with k < n are generated by the \n(u/v)— enhanced 
terms in the right-hand side of (]5.27jl . 



42 



Let us demonstrate that the contribution of scalars to the two-loop mixing matrix 7^ fc is 
fully consistent with the conformal symmetry constraints (j6.30j) and (J6.32)) . One notices that the 
scalars contribute to the right-hand side of ()6.32|) through the one-loop anomalous dimension, 
Eq. (J6.8j) . and the one- loop beta-function 

\n s N c + ... , & = -- 
2 6 

where n s = 2(jV — 1) and ellipsis denote the remaining, n s — independent terms. Substituting 
these expressions into ()6.32j) . we can predict the contribution of scalars to the off-diagonal part 
of the two-loop mixing matrix 



7i 0) = ^sN c + ..., fa = --n s N c + ... , (6.40) 



(l)sc 
7nfc 



n>k 



7f-7l 0) )^(-y + y)^ c , (6.41) 



where is given by (jfj.Hjl and the coefficients d n k were defined in (16.33)1 . A straightforward 
calculation of the integral in the left-hand side of 1)6.39)1 leads to the same expression. 

6.1.5 Anomalous dimension of twist-two operators 

Let us finally turn to the calculation of the eigenvalues of the two-loop mixing matrices in QCD 
and SYM theories. We remind that the mixing matrix has a triangular form 1)6.24)1 and, therefore, 
its eigenvalues coincide with the diagonal matrix elements, 7n(g 2 ) = 7nn(<? 2 )- The off-diagonal 
matrix elements correspond to the admixture of operators with total derivatives and they do not 
affect 7„(g 2 ). One can easily get rid of this redundant mixing by neglecting the operators with 
total derivatives. This can be rigorously done by going over from the operators to their forward 
matrix elements over some reference momentum state, the so-called forward limit 

0nn(O)->(p|0„„(O)b>. (6.42) 

In distinction with O nn (0) its forward matrix elements have an autonomous scale dependence 
since the operators with total derivatives do not contribute (p\O nn (0)\p) . In the momentum 
representation, the forward limit amounts to setting the total momentum equal to zero, J2k Uk = 

To two-loop order, the eigenvalue of the mixing matrix can be written as 

7„( S 2 ) = |U°> +(^) 2 ^ ) + ---- < 6 - 43 > 

Here is the same in the MS— and DR— schemes while is scheme- dependent. In QCD, the 
eigenvalue of the two-loop evolution kernel ()5.19)) can be extracted from the available expression 
for the anomalous dimension of the quark transversity operator in the MS scheme jSZ| 

7 <°) = C F {±S 1 (n)-Z) (6.44) 

7 «.« = C 2 j_3 + + 16 £_ 2ji(n) _ 8S3(n) _ 85 _ 3(n) _ 8Sl {n)(S 2 {n) + 25_ 2 (n)) 

f 17 134 22 

+ C F N C |- - + — S 1 (n) - -S 2 (n) - 8SL 2>1 (n) + AS 3 (n) + AS. 3 (n) + 8S 1 (n)S^(n) 

f 1 20 4 1 

+ C F n f |- - -S 1 (n) + -S 2 (n)j , (6.45) 
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for odd n. It is expressed in terms of the harmonic sums 

n+1 C_i_i \k n+1 C_i_i Nfc 

= E V- ' ^ m(n) = E - 1} • (6 - 46) 

fc=i /c=i 

The one- loop result (|fi.44j) coincides with ()6.7j) . The two-loop result can be reproduced by 
calculating the mixing matrix in the basis of conformal operators, 7^^, and, then, finding its 
eigenvalues. 

To obtain the eigenvalues of the two-loop mixing matrix in SYM case, we have to "supersym- 
metrize" the QCD expression ()6.44j) and (J6.45J) by adjusting the values of the color factors, then 
adding the contribution of scalars present in SYM theories with extended supersymmetry and, 
finally, transforming the entire result to the DRED scheme. The first of these steps is trivially 
accomplished via the identification 

C F ->N C , n f ~^N c Af. (6.47) 

The anomalous dimension obtained via this procedure will be denoted ^ CD > MS . We remind that 
in SYM with M supercharges, the gaugino fields belong to the adjoint representation and their 
total number is related to the number of scalars as rif = 1 + n s /2 = M . So we set it to this value 
in what follows. 

Let us now add the contribution of scalars. To one- and two-loop orders it comes from 
diagrams shown in Fig. [T] (c) and Figs. H3 El respectively. These diagrams have been calculated 
in the previous section as a part of the evolution kernel for three-gaugino operators. For the 
twist-two operators, the only difference is that the pair of gaugino fields carries zero color charge 
which affects the overall color factors. Using our findings from the previous section, we get for 
the contribution of scalars in the MS scheme 

Combining this expression together with the QCD result "supersymmetrized" via (|6.47|) . we get 
the anomalous dimension of the twist-two gaugino operator (|2.17j) in the MS— scheme 

7r (g 2 ) = in QCD > m (g 2 ) + in' m (g 2 ) ■ (6.49) 

As a final step we have to transform this result to the supersymmetry preserving dimensional 
reduction scheme. This scheme transformation was described in detail in Sect. 5.3.1 for a general 
L— particle operator. Applying the formulas ()5.53j) for L = 2 we obtain the anomalous dimension 
in the dimensional reduction scheme as 

iT(gk) = 7f (gk- f fp) - ^ (fp) 2 + otf) . (6.50) 

Combining together ()6.49|) and ()6.50|) . we find the one- and two-loop corrections to the anomalous 
dimension ()6.43|) in the SYM theory as functions of the number of supercharges M 

7 i 0) =N c {U-A + AS l {n)}, 

7 (i)^ = N 2 C {-(M - 4)(jV - 2) - 4 [A/" - 4 + 2S 2 (n) + 2S- 2 (n)] S^n) 

-AS 3 (n) - AS- 3 (n) + 85_ 2]1 (n)} . (6.51) 
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Here, the Af— dependence accounts for contribution from n s = 2{Af — 1) scalar and rif — Af 
gaugino fields. The one- loop result matches ()6.8|) . 
The following comments are in order. 

In SYM theory, the maximal helicity twist-two operator ()2.17|) is a component of a supermul- 
tiplet of twist-two operators. Supersymmetry implies that all components of the supermultiplet 
have the same anomalous dimension, Eqs. ()6.43|) and (j6.51|) . In the maximally supersymmetric 
Af = 4 SYM theory, all twist-two operators belong to a single supermultiplet and, as a con- 
sequence, their anomalous dimension is given by the so-called universal anomalous dimension. 8 
One can verify that for Af = 4 our result for the anomalous dimension, Eqs. ()6.43|) and (j6.51|) . 
agrees with the one of Refs. [5T| I52j. 

The expression for the two- loop anomalous dimension, Eqs. (|6.43|) and (j6.51|) . can be ana- 
lytically continued from nonnegative integer n to the entire complex n— plane. The value of the 
anomalous dimension for n = —1 is of a special interest. Formally, it corresponds to an anomalous 
dimension of a nonlocal operator X^ A d^X B ^ with canonical dimension 2. However, supersym- 
metry implies that the same quantity determines the anomalous dimension of a local operator, 
given by the product of two complex holomorphic scalar fields. The latter operator is protected 
by supersymmetry and its anomalous dimension is given to all orders by the beta-function of the 
corresponding SYM model jHSJ- This leads to the following relation 

72-i(s§h) = /W<7rk) ■ ( 6 - 52 ) 

Indeed, one verifies using (J6.46j) that for n — — 1 all harmonic sums vanish in the left-hand side 
of (j6.51j) and the relation ()6.52|) holds true to two loops. 

We argued in Sect. 15.4.31 that the AT— dependence of the evolution kernel for three-gaugino 
operators can be factored out into c-valued normalization factors. One can verify that the same 
relation (j5.64j) fulfilled by the anomalous dimension of the twist-two operators, Eqs. (j6.43|) and 
(J6.51|) . Moreover, taking into account ()6.52|) one can write this relation in a concise form as 

[ln(g 2 )-P SYM (g*)] JV = (l-(^-l)^) bn(g 2 )-P SYM (9 2 )} M=1 + 0(g e ), (6.53) 

where we introduced the subscript to indicate the value of Af = 0, 1, 2, 4. We conclude that the 
anomalous dimensions of twist-two operators f)2. 17j) in all SYM theories are related to each other 
to two- loop accuracy as in (j6.53|) . Taking into account that j3 SYM (g 2 ) = for Af = 4, the same 
relation can be written as 

[7n(<7 2 ) - (3 SY m(9 2 )] m =(l-(Af- 4)^) [7nG? 2 )V=4 + 0(g e ) . (6.54) 

Let us examine the asymptotics of the anomalous dimension at large n. For n ^> 1 the leading 
contribution to (|6.51|) comes from terms involving S\(n) ~ Inn. In this way one finds that the 
anomalous dimension scales logarithmically at large n 



ln{g 2 ) = 4 



g 2 N c fg 2 N c " ' 



?tt 2 V 8vr 2 



O 



lnn + e>(n°). (6.55) 



8 In SYM theory with Af < 4 supercharges there are two different supermultiplets of twist-two operators and, 
as a consequence, there are two universal anomalous dimensions |24) . The maximal helicity operators belong to 
one of the supermultiplets. 
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Comparing this relation with ()5.57|) one observes that the prefactor coincides with the cusp 
anomalous dimension leading to 

ln {g 2 ) = 2r cusp (A) In n + 0{n°) , (6.56) 

with the coupling A = g 2 N c / (8n 2 ) defined in (|4.1(J|) . One can show that this relation holds true to 
all loops with r cusp (A) modified by higher order corrections jH]. Substituting ([6.56]) into ([6.53]) 
one finds that as far as the M— dependence is concerned the cusp anomaly satisfies the same 
relation, Eq. (|6.53jl . 

6.2 Three-particle operators 

Let us now turn to diagonalization of the evolution kernel for the three-fermion operators in 
QCD and SYM theories. As for two particle operators, we start in Sect. 16.2.11 with finding 
the eigenspectrum of the dilatation operator to one-loop order. This allows us to elucidate in 
Sect. 16.2.21 the symmetries of the all-order dilatation operator and their manifestation in the 
eigenspectrum. Then in Sect. 16. 2. HI we numerically study the eigenspectrum of the dilatation 
operator ()5.46|) at two loops and consider afterwards the limits of large spin and (3q — > oo. 

6.2.1 Eigenspectrum to one-loop 

The one- loop dilatation operator has been determined in Sect. 14.21 We demonstrated there 
that it is given in QCD and SYM theories by the same universal expression. In the coordinate 
representation, the eigenproblem for the one-loop dilatation operator reads 

H (0) ^ 9 (z l5 Z2, z 3 ) = ^ q (z h Z2, z 3 ) (6.57) 

with the kernel H (0) (|4.7j) built from the pairwise interactions (j4.11j) and q enumerating the 
solutions. In the momentum representation, the same relation looks like 

J [du} 3 Y m (u\v)P q (u u u 2 ,u 3 ) = -TfP^i,^,^), (6.58) 

with the evolution kernel Y (0) (u\v) given by ([4.13]) and the integration measure [du} 3 defined in 
Eq. ([3.19]) . The advantage of the coordinate representation is that the conformal symmetry is 
manifest (see Eq. (J6.21)) ). On the other hand, the momentum representation is more convenient 
for actual evaluation of the eigenvalues 7^ 0) since it allows one to eliminate the contribution of 
operators with total derivatives by going over to the forward limit Uk = Vk = (see 
Appendix |XJ) . The eigenstates in the two representations, ^ q {zi) and P q {vi) are homogeneous 
polynomials in light-cone coordinates and momentum, respectively. They are related to each 
other through the orthogonality condition (J3.26|) . We remind that P q {vj) determines the explicit 
form of local Wilson operator ([3.27]) while ^ q {zj) defines the coefficient function accompanying 
this operator in the operator product expansion (|3.4|) . 

The dilatation operator HI , Eq. ([4.7]) . is invariant under the cyclic permutation of particles, 
P, and the interchange of any their pair, say P12, 

[P, ETj = [P 12 , H] = , (6.59) 
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with the corresponding discrete operators defined as 

Ptf,(zi, 22,33) = V q (z 2 , Z3, z x ) , Pl2*,(2l, z 2 ,z 3 ) = V q {z2,z x ,zz). (6.60) 

Since these two operators do not commute, [P, P12] 7^ 0, the eigenfunctions ty q (zi) can diagonalize 
either P or P12, but not both simultaneously. If ^ q (zi) diagonalizes P, then the eigenstates of 
(J6.57)) can be classified according to their quasimomentum 9 q , 

Fy q (z 1 ,z 2 ,z 3 ) = e i9 '>y q (z 1 ,z 2 ,z 3 ), (6.61) 

which takes three different values satisfying e 3ldq = 1 

9 q = 0,±—. (6.62) 

Then, the eigenfunctions of ()6.57|) with definite parity are given by 

«J (*) = ^(1 ± fn)^ q (z) , Pia^J («) = ±*J(z) • (6.63) 

Invariance of the dilatation operator under the discrete transformations (|6.59j) immediately im- 
plies that its spectrum is double degenerate except for the eigenstates with zero quasimomentum 
6 q = 3J. To see this one applies the eigenstate ^ q (z±, z 2 , z 3 ) to both sides of the relation 
PP12P = P12 which follows from the definition ()6.60|) . In this way, one obtains 

P (Pi2*,(*i)) = e"** Pi 2 *,(*) , (6.64) 

and, therefore, fi 2 ^ q (zi) defines yet another eigenstate with the same "energy" 7^ 0) unless 9 q = 0, 
leading to 

pW-7,(A)]*f(«j) = 0. (6.65) 

As we will see in a moment, integrability extends this property to the eigenstates with 8 q = 0. 
We would like to stress that the degeneracy of the eigenstates with 9 ? / follows from the 
symmetry of the dilatation operator under discrete transformations (J6.6Uj) and, therefore, it 
holds to all loops. 

The additional restrictions on the eigenspectrum of the dilatation operator are imposed by 
the properties of the nonlocal light-cone operators ()2.19|) and (|2.2U|) . Since both operators are 
invariant under cyclic permutations of fields, the same property should hold true for the polyno- 
mials ^(zi) entering ([3.4)1 . Therefore, their quasimomentum ought to be equal to zero, 9 q = 0. 
One can "revive" the remaining eigenstates with 9 q 7^ by assigning an additional "flavor" index 
to fermions. In QCD this procedure has a direct physical meaning (there are six different species 
of quarks) while in the SYM case it would violate supersymmetry by breaking a balance between 
bosonic and fermionic degrees of freedom. 

The one-loop dilatation operator H (0) is invariant under the conformal transformations, [L" + 
L% + L3,EI (0) ] = 0, with the SL(2) generators defined in (|6.12j) . This allows one to classify 
its eigenstates ^ q (zi) according to irreducible components in the tensor product of three SL{2) 
modules, Vj <8> Vj <S> Vj. Applying the relation ()6.14|) twice, one finds that ^ q (zi) carries the total 
three-particle conformal spin equal to J\ 23 = 3j + N with j = 1 and integer N > 0. This is the 
eigenvalue of the three-particle quadratic Casimir 

9 f = L\ 2 + Lj 3 + L\ x - 3j(j - 1) , (6.66) 
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such that 

gf = J 123 ( J123 - , (6.67) 

with J 123 = 3 + iV. On the other hand, the two-particle conformal spin in any subchannel, say 
(12), is not fixed anymore 2j < Ji 2 < 2j + N. This means that, in distinction with the L = 2 
particle operators discussed in the previous section (see Eq. (|6.19jl ). the conformal invariance 
alone does not permit us to determine the eigenstates of the three-particle dilatation operator. 
Still, it allows us to write a general expression for the eigenstates *& q (zi) as a sum over the states 
possessing definite conformal spins Ji 23 = 3j + N and J i2 = 2j + n and reduce the eigenproblem 
(16.57(1 to finding the corresponding expansion coefficients. 

The conformal symmetry also dictates that the eigenstates ty q (zi) have to be orthogonal to 
each other with respect to the SL{2) scalar product (for a review see Refs. jlHl E>4J)- in the 
momentum representation, the same condition looks like 

(y q \Vq,)= [ [dw] 3 MiM 2 « 3 P q (Mi)P q 'K) = 8 qq r . (6.68) 
Jo 

This relation is a counter-part of the orthogonality condition for the Gegenbauer polynomials 
defining the conformal polynomials for L = 2 operators, Eq. ((6.3)1 . For the eigenstates of definite 
parity, the orthogonality condition looks like 

(*J|*$> = <W. (*J|¥J> = 0. (6.69) 

The dilatation operator is a self-adjoint operator with respect to the scalar product ((6.68(1 and, 
as a consequence, its eigenvalues j q are real. 

The eigenproblem for the one-loop dilatation operator ((6.57(1 can be solved exactly thanks to 
the existence of the hidden, SL{2) invariant conserved charge 

p<°>, qf ] ] = [Lf + L« + LI gf ] = , (6.70) 

where we introduced a superscript to indicate that this charge could acquire higher order per- 
turbative corrections in coupling constant. The charge q% is defined as 

?f = ^L,£y = ^L^ 2 ], (6.7i) 

with the two-particle Casimir operators L 2 - k given by ()6.18(1 . As a result, the Schrodinger equation 
((6.57(1 turns out to be completely integrable and the one-loop dilatation operator rf ) can be 
mapped into the Hamiltonian of a completely integrable Heisenberg magnet of length L = 3 and 
spins being the generators of the conformal SL{2) group (for a review see Ref. [0]). 

According to ((6.70(1 . the eigenstates of the dilatation operator can be chosen to diagonalize 
the conserved charge 

qf % q (zi, 22, z 3 ) = q 22, z 3 ) . (6.72) 

Together with the total conformal spin J 123 = 3j + N, the charge q define the total set of 
quantum numbers parameterizing the solutions to the spectral problem (j6.57() . q = (N,q). A 
distinguished feature of q% is that it is invariant under cyclic permutations and changes the sign 
under permutation of any pair of fields 

[P,gf] = {P 12 ,gf } = 0. (6.73) 
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As far as the dependence of rf ) on the conserved charge is concerned, one finds using ()6.59j) 



H(°)(gf ) = P 12 H(°)(gf )P 12 = H(°)(P 12 gf P 12 ) = H^-gf ) , (6.74) 
that is, the eigenstates carrying opposite values of the conserved charge have the same energy 

= P12 , 7? } = 7^ , (6.75) 

where — q = (N, —q). For q = one finds from ()6.75|) that \l/ 9 =o(^) has the parity +1 and, 
therefore, \I>~ (^) = 0. For q ^ 0, it follows from ()6.72j) and (j6.75j) that the operator qf^ maps 
into each other "degenerate" eigenstates with positive and negative parity 

gf *J(zi,*2,*3) =q^(z u z 2 ,z 3 ). (6.76) 

We already observed that the symmetry of the dilatation operator under permutations of fields 
leads to the degeneracy of eigenstates with nonvanishing quasimomentum. According to ()6.75|) . 
integrability extends this property to the eigenstates with 6 q = and q ^ 0. The eigenstate with 
q = is the only one which is not paired. It has the quasimomentum 8 q=Q = and only appears 
for even N. 

The exact solution to the eigenproblem is given by the Bethe Ansatz technique. For given 
N, the eigenspectrum is determined by a set of N real numbers Ai, . . . , A^, the so-called Bethe 
roots, satisfying the system of transcendental equations 

^Y= fl F 1 ^' (n = l,...,JV). (6.77) 

k=l,k=£n 

The corresponding exact values of the "energy" , quasimomentum and the conserved charge are 
given, respectively, by 

N 9 N \ - ■ N f ■ \ 

fe=i k k=i K k=i v kj 

We remind that cyclic symmetry of the three-particle states imposes the selection rule 8 q = 0. 
The results of the calculation of the one-loop anomalous dimension 7^ 0) for < < 20 is shown 
in Fig. El Although 7™ can not be written in a closed form, one can work out asymptotic 
expressions which approximate the exact result with a good accuracy j^Sl IH El Ej 



6.2.2 Criterium of higher loop integrability 

Integrability played a vital role in the diagonalization of the one-loop dilatation operator. It 
appeared atop of the conformal symmetry and simplified significantly the solution of the spectral 
problem ()6.57|) but its origin remained unclear. The question arises whether this symmetry will 
survive to higher loops in QCD and SYM theories. If integrability were ultimately tied to the 
conformal symmetry then one should expect, based on the discussion in Sect. 16.1.31 that both 
symmetries are broken starting from two loops in gauge theory with /3(g 2 ) 7^ 0. As we will argue 
below this does not happen and integrability survives the conformal symmetry breaking. 

To two-loop order, the dilatation operator takes the form EI = Arf ) + A 2 !!* 1 ) + 0(X 3 ), 
Eq. ()3.2j) with the two-loop correction given in the momentum representation by ()5.46|) . The 
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eigenspectrum of tf ) was determined in the previous section. Thinking about EI as a quantum 
mechanical Hamiltonian, one can apply the perturbation theory to calculate the corrections to the 
eigenspectrum of the two- loop dilatation operator induced by the "perturbation" M.^\ In spite 
of the fact that H does not respect the conformal symmetry, the eigenf unctions of the one-loop 
dilatation operator provide a complete set of functions endowed with the scalar product ()6.68j) . 
In other words, the eigenstates of the all-loop dilatation operator HI can be decomposed over the 
leading order eigenstates ty q (zi). Since HI has the same symmetry under cyclic permutations of 
particles as rf ), this decomposition involves the states with the same quasimomentum 9 q . As 
before, the eigenstates of HI with 9 q ^ are double-degenerate while for q = the degeneracy 
can be lifted in general. 

Applying the conventional perturbation theory we find that the eigenstate of HI is given by 

= *«(*) + A £ ^^^ M*) + 0{\ 2 ) . (6.79) 

However, due to degeneracy of the spectrum ([6.75)1 . the expansion is not well-defined for q' = —q. 
We recall that the paired eigenstates have opposite parity, ty^fe), Eq. ([6.63)1 . The perturbation 
jjW preserves the parity, [HI* 1 ), P 12 ] = 0, and, therefore, the eigenstates with positive and negative 
parity get decoupled from each other 

(#J|H (1) |^) =0. (6.80) 

Then, in the sector with a definite parity, the energy levels of rf ) are not degenerate and the 
conventional formulas are at work. In particular, the two-loop correction to the eigenvalue of the 
dilatation operator with positive /negative parity is given by 

7 (i)^=(*J|H(i)|*J>, (6.81) 

with the normalization condition ([6.69)1 . Substituting *ffg(zi) by their expressions ([6.6 3)) the 
matrix element can be evaluated as 

(#J|H«|*±) = (^ q |H«|* q ) ± (* q |H«|*_ q ) , (6.82) 

where we took into account that (*_,|HW |*_,) = (# q |P 12 H[Wp 12 |# g ) = (^ q \M^ |*,). Thus, 
the two- loop correction to the dilatation operator lifts the leading order degeneracy, 7 q °' ) ' ± = 7 q °\ 



Eq. (EH5J), 



7 (D 1+ _ 7 (D,- 



2(^ q |e (1) |^_ q ) . (6.83) 

We expect that 7q 1 ' ) ' + = 7 q 1 ' ) ' for the eigenstates with 6 q ^ 0. Indeed, examining the identity 
(\I/q|[P, HlWjI^-q) = it is easy to see that the matrix element in the right-hand side of ([6.83)1 
vanishes for 9 q ^ 0. Thus, the two- loop correction to the dilatation operator could only affect 
the pairing of the eigenstates with the quasimomentum 9 q = 0. 

The necessary and sufficient condition for pairing to persist to two loops is (for q ^ 0) 

(^ 9 |H (1) |^_ q ) = 0. (6.84) 

It turns out that this relation implies the existence of the conserved charge 

qs = ?f + A^ 1} + 0{\ 2 ) , (6.85) 
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which has the same discrete symmetry as the leading order charge, [P, q^] = {P12, 93} = and 
commutes with the dilatation operator [g 3 , HI] = 0. When expanded in powers of A, the last 
relation reads 

[qi 0) ,B^] + [4 1 \B^] = 0. (6.86) 
Going over to the matrix elements one obtains 

(q-q')(* q \MW\* q ,) = ( 7 (°> -T^K^I^I**')- (6-87) 

One puts q' = —q, applies (j6.75j) and arrives at (j6.84j) . Eq. (J6.87)) can be used to determine the 
matrix elements (ty q \q^\^ q >) for q 7^ ±g'. For q = ±g' both sides of the relation ()6.87j) vanish 
simultaneously. The reason for this is that ()6.86|) defines the charge q% up to an addition of 
an operator commuting with the one-loop dilatation operator. This ambiguity can be fixed by 
imposing additional conditions on q% . For instance, we will show below that can be defined 
in such a way that the two-loop charge (J6.85)) has the same eigenvalues as the one-loop charge 
while its eigenstates are modified by ~ A corrections. 

Let us examine expressions for the eigenstates of two-loop dilatation operator 

= *J (*) + A \o)_ (o) = (* + AZ ) > ( 6 - 88 ) 

q'^q 7<J lq> 

where in distinction with 1)6.79)1 the sum in the right-hand side is well-defined for all q' . Here the 
notation was introduced for the operator Z 

Z|^) = E AO) (0) !*?>• ( 6 - 89 ) 



q>^ q % ~ 7, 



q' 



It is easy to verify that this operator is antihermitian, (Z)' = — Z. 

The operator (1 + AZ) has a simple physical meaning — it rotates the one- loop eigenstates 
into the two-loop ones. Let us perform the following unitary transformation of the two-loop 
Hamiltonian 

H = (1 - AZ)H(1 + AZ) = H (0) + A + [H (0) , Z]) + £>(A 2 ) . (6.90) 

By construction, EI has the same eigenvalues as the two-loop Hamiltonian H whereas the cor- 
responding eigenstates coincide with the one-loop eigenstates H|\l/q) = 7,1^,) with 7^ = 

7q°' ) + A7q 1 ' ) ' ± + 0(X 2 ). Combining this relation together with 1)6.75)1 one obtains 

(Hgf - gf i)|¥J) = g( 7 f - 7 q T )l*J> = [Mil^ ~ 7^' T ) + 0(X 2 )} |*J> . (6.91) 

If the degeneracy of eigenvalues survives to two loops, g(7q 1 ^' ± — 7q 1 ' ) ' =F ) = 0, then the operators 
EI and gg commute leading to 

[H,gf ] = [H, (1 + AZ)gf (1 - AZ)] = + C(A 2 ) . (6.92) 

We deduce from this relation that the pairing of eigenvalues allows us to construct the operator 
satisfying f£5H) and (J6~85j) 

q, = qf -\[qf\n + 0{\ 2 ). (6.93) 
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Taking into account ()6.89|) and ()6.76|) . it is straightforward to verify that the operator = 
— [qf\z] satisfies ()6.87|) and has the same discrete symmetries, Eq. ()6.73|) . as the leading term 
By construction, the operator q% commutes (up to higher order corrections) with the two- 
loop Hamiltonian. Since q% is obtained from by the same unitary transformation as in ()6.90|) . 
its eigenvalues do not receive perturbative corrections while the eigenstates are given by ()6.88|) 



We recall that the definition of is ambiguous since one can fulfill ()6.86j) by adding to g 3 an 
arbitrary function of q% . Eq. ()6.94j) corresponds to a particular choice of the 'gauge' in which 
eigenvalues of qs do not receive radiative corrections to two loops. 

In an analogous manner one finds the higher order analogue q2 of the quadratic confor- 
mal Casimir g 2 , Eq. (j6.66|) . It is determined by the same rotation matrix Z such that q 2 = 
q% ~~ -M?2 >2] + 0(X 2 ). Notice however that the charge g 2 , having an opposite permutation 
symmetry compared to q%, does not require the pairing of eigenvalues with opposite parity for its 
commutativity with the Hamiltonian EI to hold. As a consequence of our construction, we found 
a unitary transformation which gives the two-loop Hamiltonian HI which can be diagonalized 
simultaneously with mutually commuting leading order charges g 2 ° an d ■ 

We would like to stress that existence of q$ is tied solely to degeneracy of the eigenvalues 
of the two-loop Hamiltonian HI. The same relation works in the opposite direction. If the 
Hamiltonian HI possesses the conserved charge q% which is odd under permutation of particles 
then the degeneracy of the eigenspectrum follows in the same manner as at the leading order, 
Eq. (jnHJ). 

The commutativity of the rotated two-loop Hamiltonian HI with the quadratic Casimir q^ 
does not imply the restoration of the conformal symmetry. The persistence of the latter would 
require that the local Wilson operators O q (0), Eq. ()3.4j) . corresponding to the eigenstates ()6.88|) 
have to satisfy the relations similar to ijOfll . that is [£> 9 (0),L+] = and [C q (0),L°] ~ O q (0). 
Notice that HI defines radiative corrections to L° and the second relation is indeed satisfied in 
the basis (|6.88|) . The question thus boils down to the study of L + . A simple consideration 
demonstrates that the higher order corrections in theories with nonvanishing (3— function lead 
to [Cq(0),L + ] ~ P(g 2 ) (see, e.g., [HUES)- For the sake of the argument let us discuss the 
well-studied twist-two sector, which was explored in detail in Sect. 16.1.31 As we have shown 
there, for (3 = 0, one can bring the Ward identities for both L°, Eq. (I6.27|) . and L + , Eq. (I6.28j) . 
to the needed form (|6.38|) by going to the CS scheme with the transformation matrix £>, ()6.36j) 
and ()6.37|) . respectively. Thus, the rotated operator O nn verifies ()6.38|) and defines the lowest 
weight vector of the irreducible representation of the conformal group. However, in case when 
the /3— function is nonzero, while the L° operator can still be diagonalized with the rotation 
matrix B, the commutator of the conformal boost L + with the rotated Wilson operator O nn 
is proportional to (3{g 2 ) and, therefore, O nn does not represent the lowest weight. The same 
phenomenon re-occurs for three-particle operators - while one is able to diagonalize L°, with Z, 
it is not the case for L + unless (3 = 0. 

6.2.3 Numerical diagonalization of the two-loop dilatation operator 

We argued in the previous section that the degeneracy of the eigenvalues of the two-loop dilatation 
operator could serve as a criterium for two-loop integrability. Let us now use the expressions for 




(6.94) 
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the two-loop dilatation operators obtained in Sect. El to evaluate their eigenspectrum and verify 
whether the pairing of the eigenvalues present at one-loop in QCD and SYM theories survives to 
two loops. 

We recall that the two-loop evolution kernels have been evaluated in the momentum repre- 
sentation. In this representation the eigenvalue problem looks like 



where the integration measure was defined in Eq. ([3.19)1 . while the integration region over the 
u— variables is determined by the support properties of the evolution kernel Y(u\v). Notice that 
the eigenvalues 7 q (A) do not depend on the "external" momenta v = (t>i, v 2 , U3). One can explore 
this fact to simplify the actual calculation of the anomalous dimensions. Indeed, the calculation 
of the two-loop evolution kernel has been performed in Sect. Elm the kinematics < v\, v 2 , v 3 < 1 
and V\ + v 2 + v 3 = 1 . In this case, the phase space for integrated u— momenta is the one shown 
in Fig. El 

As explained in Sec. 16.1.31 for two-particle operators, the conformal symmetry is broken be- 
yond leading order in gauge coupling and so the eigenfunctions, depending on the quantum 
number N, of the dilatation operator are changed by 0(A) corrections. This symmetry breaking 
is induced by the renormalization prescription of the composite operators. In fact the breakdown 
of special conformal symmetry arises already at one- loop level and this determines via the con- 
formal constraint ()6.3(J|) the eigenfunctions of the two-loop dilatation operator. Hence, one can 
get rid of this conformal symmetry breaking at this order by means of a finite renormalization 
of the operator to one-loop accuracy. It can be shown in general that the known special confor- 
mal anomaly in the two-particle subchannels to one-loop accuracy together with the conformal 
constraint, similar to Eq. ()6.30j) . is consistent with our result (|5.46J) . 

For the (3(g) = case, the finite renormalization alluded to above yields the anomalous 
dimension matrix which is identical to the one found by means of taking the forward limit. 
Therefore, the simplest way to calculate the eigenvalues of the two-loop evolution kernel within 
the conformal subtraction scheme is to take the forward limit of the kernel, i.e., impose the 
following condition on the light-cone particle momenta 



As was already explained in Sect. 16.131 the forward limit effectively corresponds to taking forward 
matrix elements of three-particle operators. This does not change the eigenspectrum of the 
dilatation operator but drastically simplifies the numerical calculation of moments since the 
mixing with the operators containing total derivatives is automatically removed. 

In the forward limit, the evolution kernel looks differently as compared to the one in the 
off-forward kinematics Yln v n = 1- ft turns out that the two kernels are related to each other 
through the following limiting procedure [HSj. Let us set 




(6.95) 



U\ + u 2 + u 3 = f 1 + v 2 + v 3 = . 



(6.96) 





Xi + X 2 + £3 = 2/1 + 2/2 + 2/3 = r 



(6.97) 



and take the limit r 



0. Then, the kernel in the forward limit reads 



F(x\y) = LIMV(u|«) = lim — V 



xi T-xt-x 3 x 3 y 1 T-yx- y 3 y 3 



(6.98) 



r 



T 



T 



T 



T 



T 
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where x denotes now the variables x\ and 23 (with X2 = —x\ — 23) and similarly for y. Before 
the limit ()6.98j) can be performed, the support of the kernel must be extended. This extension 
procedure is described in Appendix |XJ Performing the limiting procedure ()6.98|) in both sides of 
the evolution equation ()6.95|) one obtains the eigenproblem in the forward limit in the following 
form 

J dx¥(x\y)C q ( Xl ,x 3 ) = — 7q(A) C q (yi,y 3 ) , (6.99) 

with dx = dxidx^. Here the eigenvalues 7 q (A) are the same as in (|6.95|) while the eigenfunctions 
are related as 

C q ( yi ,y 3 ) = lire 1/ P q T - y '- y \ Hi) = Pq {yu _ yi _ y 3 , y3 ) . (6.100) 

In the last relation we used the fact that the eigenstates of the evolution kernel are homogenous 
polynomials in momentum fractions of degree N > 0. Obviously, the polynomials C q (y 1,1/3) 
have the same property. They are uniquely specified by the set of coefficients in the expansion 
of C q (1/1,1/3) over the basis yiy 3 ~ k (with k = 0, . . . , N). The total number of these coefficients 
is N + 1 and it matches the total number of three-particle operators with N derivatives having 
nonvanishing matrix elements in the forward limit, i.e., free from total derivatives. In the basis 
Hills ~ k the forward evolution kernel ¥(x\y) is represented by a finite dimensional mixing matrix 

r - 

/ dxF(x\y) x»- n x n 3 = J2 A nkWy?- k y k 3 , (6.101) 

whose entries depend both on iV and the gauge coupling A. Then, to find the spectrum of the 
anomalous dimensions it suffices to solve the characteristic equation for the mixing matrix for 
each value of N > 

det(A- 7q (A)) = 0. (6.102) 

Evaluation of the two-loop mixing matrix in the forward limit is a straightforward but tedious 
exercise. For given N, the mixing matrix has dimension N + 1 and, therefore, the characteristic 
equation has N + 1 solutions for j q . All eigenvalues of the mixing matrix are real. This property 
is a consequence of the conformal invariance of the dilatation operator to one-loop order. In 
addition, some of the eigenstates of the mixing matrix are degenerate. To one-loop order all 
eigenstates except the one with q = are paired and the eigenspectrum contains [1 + N/2] 
different energy levels. To two-loop order, the pairing can be only lifted for the eigenstates with 
vanishing quasimomentum. We remind that the cyclic symmetry of Wilson operators selects the 
latter eigenstates. 

To begin with, one performs numerical diagonalization of the one-loop mixing matrix in QCD 
and SYM theories (see Fig. El (a))- In both cases, the one- loop dilatation operator involves the 
same two-particle kernel and, therefore, the one-loop mixing matrices are related to each other. 
One verifies that, in agreement with the exact results based on the Algebraic Bethe Ansatz, 
Eqs. (16.78)) and (|6.77|) . the one- loop eigenvalues coincide with the energy levels of the Heisenberg 
spin magnet. For given total conformal spin Ji 23 = 3j + N, there are N + 1 eigenstates. For 
odd, all energy levels are double degenerate while for N even this is true for all levels except 
the lowest lying one. The latter eigenstate carries the conserved charge q = and carries 
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Figure 9: The spectrum at one-loop order (a) and two-loop order (b) in QCD in the MS scheme. 
The eigenvalues with 9 q ^ are removed from the spectra for better comparison with SYM cases 
where only zero quasimomentum states are physical. 



the vanishing quasimomentum 9 q = 0. For eigenstates with quasimomentum 9 q = ±§7r the 
degeneracy is a consequence of the symmetry of the evolution kernel under discrete symmetries. 
The double degeneracy of eigenvalues with zero quasimomentum implies integrability of the one- 
loop Schrodinger equation ([6.57)1 . For a given N, the total number m of the eigenstates with 
9 q = is 

1 2 

m = - (N - 1) + -mod(AT - 1, 3) . (6.103) 
3 3 

Among them, as we said above, there is one eigenstate with q = for even N = 2k with 
k = 0,1,2, ... , while the rest are paired. 

To two-loop order, the dilatation operators in QCD and SYM theories are given by two 
different expressions and, therefore, the corresponding two-loop mixing matrices are not related 
to each other in a simple manner. Diagonalizing the latter in the forward limit, we observe the 
following remarkable phenomenon (see Fig. |U] (b) and Figs. ITU1 (b.c.d)): 

• In QCD, for the baryon operator Of un (zi, z 2 , z 3 ), the pairing is lifted for the eigenvalues 
with zero quasimomentum but indeed it persists for the eigenvalues with 9 q = ±|vr. 

• In SYM theories with arbitrary number of supercharges M, the one-loop pairing of eigen- 
values carries on to two loops for the operator O a dj(^i, z 3 ); 

We recall that degeneracy of the eigenstates with nonzero quasimomentum is not related to 
integrability. 

Presenting the two-loop expressions for the eigenvalues it proves convenient to introduce 

A 7 f(A) = tJ(A) - 3T(A) , (6.104) 

where 3T(A) is the anomalous dimension of the local three-fermion operator without derivatives, 
N = 0. The anomalous dimension defined in such a way satisfies the normalization condition 
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(A) = for N = 0. We recall that the anomalous dimension T(A) has the perturbative 
expansion (j3.15|) with one- and two-loop corrections in the MS scheme defined in (j4.12|) . ()5.44|) 
and (|5.45jl . For SYM theories, we use the DR scheme with r^P 11 related to (|5.45j) as 



,(1),DR 



r (i) _ I r (o) 

1 SYM n L SYM 

6 



+ I(4-^ = 2-ijV(JV-l) 



(6.105) 



The exact spectrum of eigenvalues A-f^ (A) up to conformal spin N = 20 is displayed in Figs. M 
and They clearly demonstrate the lifting of the degeneracy in QCD for the baryon operator 
()2.20j) compared to TV-extended SYM where it persists for all eigenstates. 

The first non-trial example of the pairing phenomenon arises for N = 3 and demonstrates 
the main trend, recurring for higher N th . For N = 3 there are four different eigenstates. Two of 
them have the quasimomentum 9 q = while for the remaining two possess the quasimomentum 
equal 6 q = ±2n/3. As was already explained, the latter eigenstates are degenerate to two loops 
but the corresponding Wilson operators are not cyclically symmetric. The pair of the eigenstates 
with 9 q = is degenerate in SYM theories for arbitrary N while in QCD the degeneracy is lifted 
to two loops. 

To illustrate the phenomenon, we first present the explicit expressions for the anomalous 
dimensions for N = 3 in the M = 1 SYM theory in the DR scheme 9 
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A 7 ±(A) 
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where (3q/Cr = 6, and in QCD in the MS scheme 



A7i + (A) 
A 7 f(A) 
A71KA) 



A4 + A 



A4 + A' 
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(6.107) 
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192 C R 



with Po/Cr = y — Nf. Here A7j ± (A) and A 7 ^(A) correspond to the eigenstates with 9 q = 
and 6 q = ±27r/3, respectively. Thus, for quarks in the fundamental representation, the two 
eigenstates with 6 q = have different eigenvalues starting from two loops, A7j + ^ A7f. We 
remind that in QCD the three-quark baryon operator only exists for N c = 3, so that large- N c 
counting is unavailable and the non-planar diagram in Fig.0](b) contributes and partially leads to 
breaking of integrability in QCD at two loops. One immediately sees from the above equations 
that though the degeneracy is generally lifted in QCD, the (3q— terms of the would-be-paired 
eigenvalues do coincide. In SYM they do not break integrability either. The fact that (3q— terms 
preserve degeneracy of the eigenstates implies that the two-loop integrability is not related to 
conformal symmetry. We shall return to this issue in Sect. 16.2.51 



9 In the MS scheme these anomalous dimensions coincide with Eq. (17) of Ref. ^Hj- Though, the anomalous 
dimensions A7 : ^ are not "physical" for cyclically symmetric operators, we present them to demonstrate that the 
degeneracy is always present for eigenvalues with nonzero quasimomentum. 
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Figure 10: The spectrum at one-loop order (a) and two-loop order in SYM theories in the DR 
scheme for Af = 1 (b), Af = 2 (c) and Af = 4 (d). Only eigenvalues with zero quasimomentum 
are displayed. 
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As a next step, we examine the expression for the anomalous dimension of paired eigenstates 
with vanishing quasimomentum, Eq. (J6. 106)1 . in SYM theory with M supercharges in the DR 
scheme 

,(A) = A4 + A 2 r ' )j 



12 



-AN 



(6.108) 



Starting from two loops it depends linearly on the number of supercharges M and satisfies ([5.64)1 . 
It is easy to see that the two loop correction to A7 I ± (A) is positive for J\f = 1 and negative for 
M = 2 and M = 4. As can be seen from Fig. El the same is true for higher values of conformal 
spin N. 

The two-loop correction to (J6. 108)1 receives contributions from both two- and three-particle 
kernels, Eq. ([5.46)1 . and the question arises whether each of these contributions preserves inte- 
grability separately. To this end, one multiplies three-particle kernels in ()5.46j) by some factor 
£ and examines the £— dependence of the resulting expressions for the anomalous dimension. In 
this way, one finds 



A 7 +(A) = A4 + A 2 
A 7 f (A) = A4 + A 2 



18 s 36 



(6.109) 



For £ = 1 one recovers ()6. 108)1 while for £ ^ 1 the degeneracy is lifted. Thus, the two-loop 
integrability in SYM theories only arises in the sum of two- and three-particle contributions. 10 

Let us present the exact two-loop eigenspectrum for a few lowest anomalous dimensions with 
9 q = in SYM theories with M supercharges. For odd conformal spin N all eigenvalues with 
zero quasimomentum are paired and we denote the corresponding two-loop anomalous dimension 
matrix as A 7 ^. For even N there i: 
we denote as A 7 jy. Then one finds 



matrix as A 7 ^. For even N there is an additional unpaired eigenstate with q = whose energy 
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(6.110) 



Notice the change of the sign of two-loop correction to the anomalous dimensions as one 
increases the amount of supersymmetry present in the model. While for QCD and M = 1 SYM, 
the eigenvalues are strictly positive, the trend starts to change for M = 2 SYM, where the low- 
spin anomalous dimensions are negative and decrease and then for higher iV the behavior starts 
to change and the spectrum recovers to positive values at higher conformal spins. For M = 4 
SYM, the spectrum is entirely negative. This change in the tendency can be traced back to 



10 One may also wonder whether integrability would persist if one merely iterated the twist-two twodoop kernel 
H_l=2 in all pairwise subchannels. A straightforward analysis demonstrates that the answer to this question is 
negative and the degeneracy of eigenvalues with zero quasimomentum is lifted in this case. 
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the linear in M term in the right-hand side of (j6.110j) which in its turn comes from (JHJ2J). It 
provides a large negative contribution to the anomalous dimensions for M = 2 and M = 4, which 
competes with the positive one from other terms. 

6.2.4 Large-spin limit 

One can see from Figs. El and El that the anomalous dimensions of the three-particle operators 
occupy a band 

7min(A) < 7<j(A) < 7max(A) . (6.111) 

Let us demonstrate that both boundaries scale as ~ In iV at large conformal spins N ^> 1 with 
certain computable overall coefficients. 

It is well known that the asymptotics of the anomalous dimensions at large conformal spins 
is governed by the terms in the evolution kernel which are singular for Vj — > Uj. This correspond 
to the situation when the wavelength of exchanged gluons is large. The behavior of the two-loop 
evolution kernels for Vj — > Uj was scrutinized in Sects. 15.4.21 and 16.1.51 We have found in section 
15.11 that the three-particle irreducible part of the three-particle kernel is not singular in this 
limit and, therefore, the leading asymptotic behavior of the anomalous dimensions for N ^> 1 
emerges entirely from the two-particle irreducible kernels (|5.55|) . According to ()5.56|) . the all-loop 
evolution kernel H(A) is given for Vj — > Uj by a remarkably simple expression 

H(A) = ^r cusp (A)HW + ... (6.112) 

where rf ) is the one-loop evolution kernel and ellipses denote terms subleading for Vj — * Uj. 
We would like to stress that this relation captures the leading ~ In iV— terms in the anomalous 
dimension and neglects subleading ~ iV° corrections. 

The fact that H(A) reduces to the one-loop evolution kernel accompanied by a multiplicative 
factor of r cusp (A) implies that the all- loop dilatation operator is integrable in the large-spin limit 
both in QCD and SYM theories. The underlying integrable structure is the one corresponding 
to the classical Heisenberg magnet jUj with the Hamiltonian which can be cast into the form 

H(A) = T cusp (A) lng^ + 0(N°) , (6.113) 

with <23° c \ being the classical counterpart of the quantum integral of motion This leads to 
the following asymptotic behavior of the anomalous dimension 

7<j(A) = T cusp (A) \nq + 0(N°) . (6.114) 

For N ^> 1, the conserved charge scales in the lower part of the spectrum as q ~ 0(N 2 ) and in 
the upper one as q ~ 0(iV 3 ) so that 

7min (A) = 2r cusp (A)lniV, 7max (A) = 3r cusp (A) InJV . (6.115) 

These relations are in agreement with the results of Ref. 03]. So far we studied the dilatation 
operator for nonlocal light-cone operator of length L = 3. It worth mentioning that the relation 
(I6.112j) holds true for arbitrary L. 
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6.2.5 Large (3 — limit 

As follows from ([6.106)1 and (J6.107)) . the corrections to the two-loop anomalous dimensions coming 
from the terms involving the (3— function preserve integrability. This may appear counterintuitive 
if one lands on the assumption that existence of integrability is intertwined with conformal 
symmetry. However, our findings demonstrate that the two phenomena are not related to one 
another. 

We argued in Sect. 15 .412*1 that one can enhance the contribution of conformal symmetry break- 
ing terms ~ (3q by going over to the formal limit (3q — > oo with A/?o = fixed. In this limit, the 
all-loop dilatation kernel is given by one-loop Feynman diagrams with the gluon propagator 
dressed by one-loop vacuum polarization bubble insertions. The resulting expression for the 
dilatation operator, Eq. ([5.62)1 . is proportional to the one-loop dilatation operator in which 
the "bare" conformal spin of quark/gaugino fields j — 1 is replaced by its renormalized value 
j = l+/W/(327r 2 ). 

The same effect can be understood within the e— expansion as follows "J"jc"j. According to ([4.6)1 . 
for D = 4 — 2e, the beta-function receives an additional contribution ~ e. Therefore, there exists 
a value of D = D c for which the D c — dimensional /3-function vanishes and the theory becomes 
conformal for this space-time dimension. It follows from (|4.6|) that this occurs for 

ec = ^. (6.116) 

Since the canonical dimensions of fields depends on D, it will propagate into the value of the 
conformal spin, Eq. (J6.13)) . For quarks/gauginos, one gets 

,^^-^i.M, ( , 117) 

where in the last step we used the perturbative expansion ([2.7)1 . This precisely matches our 
finding in Eq. ([5.62)1 for the renormalized conformal spin. Thus, the evolution kernel in the 
large— /3q limit coincides with the evolution kernel in the conformal invariant field theory living 
in the "critical" D c — dimensions. 11 

We conclude that in the large— /So limit, the all-loop dilatation operator H(A) is given by the 
one-loop dilatation operator with the shifted conformal spin ()6.117[) . We recall that the latter 
operator coincides with the Hamiltonian of the Heisenberg SL(2\ M.) spin chain. From the point 
of view of integrable models, renormalization of the conformal spin corresponds to deformation 
of unitary irreducible spin j = 1 representation of the SL(2) ~ SO(2, 1) group to yet another 
irreducible but not unitary representation with spin j = l+/3 <? 2 / (327r 2 ). One can easily construct 
the "all loop-order" integrals of motion (72(A) and 93(A) for the corresponding system by merely 
substituting j in Eqs. ([6.66)1 and ([6.71)1 by its renormalized value ([6.117)1 

<?n(A) = q^\ 1+ , 0ff2 = e^q^e-^ (n = 2,3) , (6.118) 
where the notation was introduced for the operator 

Z = A^. (6.119) 

4c R dj 



11 



We are grateful to A. Manashov for proposing this interpretation. 
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Expanding ()6.118|) in powers of A one reproduces (|6.93)l . 

The large (3q— limit provides a nontrivial example of all- loop dilatation operator and allows one 
to trace the deformation of integrable structures to higher loops. As in the previous section, the 
above analysis can be immediately generalized to the nonlocal light-cone operators of arbitrary 
length L. 

7 Discussion and conclusions 

In the present work we have derived the two-loop dilatation operator in the noncompact SL(2) 
sector of QCD and SYM theories with M = 1,2,4 supercharges. We have chosen the three-quark 
and three-gaugino maximal-helicity operators, respectively, as a case of study. The choice was 
driven by the closure of the sector under two-loop renormalization and the simplicity of calcu- 
lations involved. Supersymmetry allows one to extend our findings to the entire supermultiplet 
embodying operators with different particle content. 

It is known that the dilatation operator in the aforementioned gauge theories is integrable to 
one-loop order in the sector under consideration. The major motivation for the current study 
was to verify whether integrability survives to two loops. To this end we had to use a guiding 
principle which will definitely point to its existence. We have used a criterion of double degeneracy 
of eigenvalues of the dilatation operator with zero quasimomentum. The eigenvalues with nonzero 
quasimomentum are always degenerate due to invariance of the dilatation operator under discrete 
symmetries — cyclic and parity permutations. We found that in QCD the degeneracy is lifted while 
in SYM theories it persists to two loop order for arbitrary number of supercharges M . In the latter 
case, the spectrum of the two-loop dilatation operator depends on Af. We demonstrated that 
for the light-cone operator of arbitrary length L = 2, 3, . . . the A/"— dependence has a remarkably 
simple form, Eq. ()5.64|) . This implies that the dilatation operators in the SL(2) sector of all 
SYM theories share the same integrability properties. 

Making use of the degeneracy to two loops, we constructed the conserved charges q 2 (A) and 
g 3 (A) which commute with the two-loop Hamiltonian and among themselves. The first term in 
the perturbative expansion of these charges, q% arid q% , coincide with the integrals of motion 
of the SL(2) Heisenberg magnet of length L = 3 and spin j = 1 HI El - As such, they are given 
by conformal invariant polynomials in two-particle SL(2) Casimirs L\ n (with k ^ n = 1,2, 3), 
Eq. (|6.18|) . and can be obtained from the auxiliary transfer matrix in a standard manner |21j . 
The interpretation of higher order corrections to 52(A) and 93(A) is an extremely nontrivial task 
both due to breaking of conformal symmetry in two loops and due to the fact that the SL(2) spin 
magnet is noncompact. Namely, the two-particle Casimirs L\ n are given by infinite-dimensional 
matrices in the conventional oscillatory basis and their powers are independent on each other. 
This should be contrasted with the compact SU(2) spin— 1/2 magnet in which case the spin 
operators are given by Pauli matrices and all powers of two-particle Casimirs L\ n are expressed 
in terms of unity operators and its first power. As we noticed in Sect. 16.2.21 another possibility 
to identify integrable structures to two loops could be to perform a unitary transformation of 
the dilatation operator H(A) and conserved charges g n (A) with the operator 1 + XL. It removes 
higher order corrections from the conserved charges and yields the Hamiltonian HI = rf ) + Atf 1 ) 
which commutes with (with n — 2, 3) 

[gf,H(A)] = + O(A 2 ). (7.1) 
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Since the one- loop Hamiltonian is a function of the conserved charges rf ) = tf°^(g^ , q% ), we 
conclude that the leading and the next-to-leading Hamiltonians are mutually commuting 

[H (0) ,H (1) ] =0. (7.2) 

Thus, is a "higher" Hamiltonian of the SL(2) Heisenberg spin chain. In principle, it can 
be constructed from by applying the boost operator [57, 58j B = Ylk ^^kk+v that is 
[B, . . . [B, HW]], and forming their linear combination with suitably adjusted coefficients. In 
order to meet the parity requirement, only even powers of B may appear in the expansion. 

We demonstrated that the integrability is not tied to the survival of the conformal symmetry 
at two-loop order. On the contrary, in QCD case, the beta-function term in the dilatation 
operator preserves the pairing of eigenvalues contrary to the conformal invariant part of the 
Hamiltonian. This phenomenon carries on to SYM models where the conformal invariant part is 
integrable as well. 

We found that the nonplanar diagrams provide vanishing contribution to the dilatation op- 
erator for L = 3-particle operators in SYM models. For L > 3, they induce next-to-nearest 
neighbor interactions at leading order, which are however suppressed by 1/N% with respect to 
the leading nearest-neighbor contributions. These corrections explicitly break integrabil- 

ity but this phenomenon can be restored by going to the large- iV c limit. Therefore, the large- N c 
limit is a necessary condition for integrability to survive. Note that to two-loop order there are 
nonplanar Feynman graphs (like in Figs. 121 (b) and0] (b)). These diagrams were not considered 
presently since they are accompanied by a vanishing color structure for L = 3 and also die out 
in the multi-color limit for L > 3. For the QCD baryon operator the multi-color counting is not 
applicable since the corresponding operator exists only for N c = 3. Even if one sets by hand 
the nonplanar two-particle graph to zero, the pairing would not be restored due an intricate 
interplay between color factors of different Feynman diagrams in the fundamental and adjoint 
representations. 

The two-loop dilatation operator has a rather complex structure in the momentum represen- 
tation. It is however in this representation that higher-loop computation machinery is developed 
the best. It would be interesting to translate the obtained expressions into the coordinate, light- 
cone position representation and, eventually, write down the two-loop dilatation operator in terms 
of light-cone superfields as it was done to one-loop order in Ref. [21|. m this form, all symmetries 
of the dilatation operator are manifest and they may well hint on its possible structure to higher 
loops. 
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A Forward limit 



In this Appendix we describe the limiting procedure which allows one to obtain the three-particle 
kernel in the forward limit u\ + u 2 + u 3 = v 1 + v 2 + v 3 — from the expression ()5.7|) defined on 
the simplices ()4.16j) and (|4.18j) . 

As the first step, one extends the expression for the evolution kernel (jf>.7)) to the entire domain 
—00 < u\ + u 2 + u 3 = v 1 + v 2 + ^3 < 00 by making use of the support properties of defining 
Feynman diagrams. In the two-particle subchannels, the procedure is well-known and it amounts 
to substituting 

e(vj-Uj) -> Q(u a ,v a ) , (A.l) 

with generalized step function defined in ()4.14|) . The analogous procedure can be formulated for 
the three-particle irreducible diagram too. As we found in Sect. 15. ll the three-particle kernel ()5.7|) 
is given by the sum of various F— functions which are accompanied by the following step- function 
structures 

tf£ 3 = 9{u x )9{u 3 -v 2 - v 3 )9(l - Ul - u 3 ) , (A.2) 

df 2 \ ee e( Ul )e( Vl - Ul )e(u 3 - v 3 )9(v 2 + v 3 - u 3 ) , (a.3) 

= _ v j 9 ( U3 _ v ^ e{ i _ Ul _ U3 ) . ( A .4) 

The remaining two structures ^[23 for i — 5,6 are found from the above by permutation of 
arguments as in Eq. (j5.8j) . It turns out that the generalized step-function structure can be 
restored through the following substitutions 

#£3 - & < 3 3 = ®(uuv l )[0(u 2 ,l-u l -v 3 )-e(l-v l -u 3 ,v 2 )} , (A.5) 

^123 - ©SS = ©K, - «1 - "3, U 2 ) , (A.6) 

<23 -> = ©K 1 - «l - «3)9(1 - ^ - W 3 , « 2 ) • (A.7) 

Having defined a generalized three-particle kernel in this way, the reduction to the forward limit 
can be easily performed with the help of (|fi.98j) . In the following we assume without loss of 
generality that v\ = yi/r and v 3 ee y 3 /r are positive in ()6.98|) . Consequently, in the forward case 
we find for the building block ([A. 1)1 : 

. 9(x a )8(y a — x a ) , a = 1,3 

UMQ(u a ,v a )={ V a) {ya ah . (A.8) 

-9(-x 1 - x 3 )9(x 1 + x 3 - yi - y 3 ) , a = 2 

The forward limit of Eqs. (|A.5J) . (|A.6|I . and (|A.7|) can be easily calculated and, one finds 

LIM eg 3 = 6{x x ) 9( yi - Xl ) [9(x 3 + yi ) - 6{ Xl + x 3 )} 9(y 3 - x 3 ) , (A.9) 
LIM eJI = -9{ Xl ) 9{y x - Xl ) 9{x 3 + Vl ) 9{y 3 - x 3 ) , 
LIM ejJs = 9{x 3 ) 9( Vl - Xl ) 9{x x ) 9(y 3 - x 3 ) . 
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B Subtraction terms 

The two-particle kernel ()5.28|) receives contribution from the subtraction terms (|5.29J) involving 
the integral of three-particle irreducible kernel 



W{u 1 ,u 2 \v 1 ,v 2 ) = / du' 3 ¥^ 3 (u x ,l - m - u' 3 ,u' 3 \vi,v 2 ,v 3 ) 



The evaluation of the integral is straightforward and leads to 

W(ui, U 2 \vi, V 2 ) = 9(ui)9(u X -V 2 - V U 2 \vi, V 2 ) 

+ 9{u 2 )0(v 2 - u 2 )s 2 {u x , u 2 \v x , v 2 ) + 9(ui)9(vi - Ui)s 3 (ui, u 2 \v x , v 2 ) 

where the notation was introduced for the functions 
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L(ui - v 3 , v 3 , v 2 ) + L(ux - V!, vx, v 3 ) + L(ui - v x ,vi, v 2 ) 



1 

2v~ 2 



4 -In 



2 V 2 - U 2 

V 2 



(B-3) 



(B-4) 



- L(u x - vi,v 2 , v 3 ) + L{u x - vi, vi, v 3 ) + L(ui - v x ,Vi, v 2 ) 



- 2f(u 2 ,v 2 ) 

s 3 (ui,u 2 \vi,v 2 ) = f{ui,Vi) 
V 2 + Vi 



V 2 + Vi 
V 2 Vi 



Li 2 



u 2 



iln 2 ^ + ln^ln^ + 2Li 2 

2 v 2 v 2 vi 



V 2 + Vi 
Vl - Ui 



Li 2 — 



+ 



V 2 Vi 



In — In h Li 2 

v 2 V 2 + Vi 



vi - Ul 
u 2 



u 2 



+ Li 2 



(B.5) 



v 2 



V 2 + Vi 



C(2) 



with the /— and L— functions defined in (j5.14j) and ()5.32|) . respectively. 

The subtraction term ()5.29|) provides a genuine contribution to the two-particle irreducible 
kernel V12, Eq. (|5.28p . Similar contribution also comes from the gluon kernel (J5.25|) which involves 
the Vgv (u, v )— term. It is instructive to examine the sum of three terms 



AVi 2 (u\v) = W(ui, u 2 \vi, v 2 ) + W(u 2 , ui\v 2 , vi) + 



Ui + u 2 



-V. 



(i) 



Vl 



Ui +U 2 Vi+ V 2 



. (B.6) 



Taking into account Eqs. (jB.2|) - (|B.5|) and ()5.24|) we find that various terms cancel against each 
other leading to 

AVi 2 (u\v) = 9(u 2 )9(u 2 - vi - v^&F^ui^^vi^^ + 6{ui)9{vi - ui)AJ z 2 (ui,u 2 \vi,v 2 ) (B.7) 

Ui <-> u 2 



+ 



Vi <-> V 2 
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with (t> 3 = 1 — v i — v 2) 

AJ 7 1 (u 1 ,u 2 \v 1 ,v 2 



u x u 2 



L(u 3 - v 1} v 3 , v 2 - u 2 ) - L(u 3 - vx,v 3 , v 2 ) (B.8) 



v 3 (v 2 - u 2 ) 
+ L{u 2 - v 2 , v 2 , v 3 ) + L(u 2 - v 2 , v 2 , v 1) , 
A^2(«i, u 2 \vi, v 2 ) = -L(u 2 - v 2 , vx,v 3 ) + L(u 2 - v 2 , v 2 , v 3 ) + L(u 2 - v 2 , v 2 , vi) 



31n^ + ln 2 ^ 



(B.9) 



We observe that it does not involve dilogarithms. 
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